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GLOSSARY  OF  ABBREVIATIONS 
AND  SYMBOLS 


CARDIN 

CFG 

CPU 

DFT 

DSRN 

ECM 

ECCM 

EDITOR 

G[  ] 
r'[  ] 

FFT 

GCOS 


GHz 

(Gc) 

"[ 

] 

HIS 

635 

IF 

LSB 

TSS  subsystem  which  allows  submission 
of  batch  jobs  from  a remote  terminal. 

Simulation  configuration 

Central  Processor  Unit.  The  central 
and  arithmetic  portion  of  a computer. 

Dirac  delta  function.  Unit  impulse  at 
time  equal  to  t^ 

Discrete  Fourier  Transform 

Data  Set  Reference  Number 

Electronic  Counter  Measures 

Electronic  Counter  Counter  Measures 

TSS  subsystem  which  facilitates 
modification  or  correction  of  TSS  files 

Fourier  transform  operator 

Inverse  Fourier  transform  operator 

Fast  Fourier  Transform  algorithm 

Honeywell  General  Comprehensive 
Operating  Supervisor 

Gigahertz 

Hilbert  transform  operator 
Honeywell/GE  635  computer  system 
Intermediate  Frequency 

Least  Significant  Bit 
A— 1 


LTI 

MCC 

MRN 

MSB 

MTI 

ns  (or  NS) 

PASS 

PRMFL 

RADSIM 

RCS 

RF 

RV 

STF 

s*(t)  or  S*(f) 
s*(t) 

s*(f) 

s(t) 

S+(f) 

S_(f) 
s (t) 


Linear  Time  Invariant 
Module  Classification  Code 
Module  Reference  Number 
Most  Significant  Bit 
Moving  Target  Indication 
nanosecond  (10-^  second) 

Pass  through  a simulation  configuration 

Permanent  disc  storage  file 

Radar  System  Simulation  Model 

Radar  Cross  Section 

Radio  (Radiated)  Frequency 

Reentry  Vehicle 

Stimulus/Transfer  Function 

Complex  conjugate  of  s(t)  or  S(f) 

A function  composed  of  samples  of  s(t) 

The  Fourier  transform  of  a sampled 

function,  such  as  s (t) 

* 

Hilbert  transform  of  s(t) 

- Positive  frequency  components  of  the 
signal  s(t) 

Negative  frequency  components  of  the 
signal  s(t) 

The  function  of  time  which  corresponds 
to  the  positive  frequency  components 
of  s(t) 


•< 

i 


¥ 
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S_(t) 

- The  function  of  time  which  corresponds 
to  the  negative  frequency  components 
of  s(t) 

Sj(t) 

In-phase  video  channel  representation 
of  the  bandlimited  signal  s(t) 

Sq(  t) 

- Quadrature  video  channel  representation 
of  the  bandlimited  signal  s(t) 

TSS 

Time-sharing  system 
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APPENDIX  B 


SOFTWARE  MODULES  USED  IN 
THE  RADAR  SIMULATION  MODEL 


The  software  modules  used  in  the  radar  simulation  module 
are  listed  herein  by  module  name  and  by  reference  number. 

B.l  SOFTWARE  MODULES  LISTED  BY  MODULE  NAME 

The  following  table  includes  all  of  the  software  modules 
listed  in  alphabetical  order  by  module  name.  Column  headings 
in  the  list  are  defined  as  follows: 

1.  Module  Ref.  Number:  This  is  the  number  by  which 

the  module  is  called  in  any  system  model  program. 

2.  Input:  Locations  where  input  data  to  module 

must  be  located  before  execution. 

3.  Output:  These  are  the  locations  where  output 

data  is  placed. 

4.  Subroutine  Name:  This  column  contains  a list  of 
subroutines  used  by  module. 

5.  Subroutine  PRMFL:  This  column  contains  a list  of 

the  permanent  file  where  the  subroutine  is 
located.  This  file  must  be  loaded  in  any  system 
model  program  that  uses  module. 

6.  System  Parameters:  The  system  parameters  are 

constants  that  are  used  in  the  module  and  are 
defined  by  module  number  301,  (SYS).  Parameters 
are  described  in  Appendix  C. 

7.  Module  Parameters:  Module  parameters  are  constants 

that  are  used  in  module  and  are  defined  in  a 
namelist  immediately  following  the  module  call  in 
the  system  model  program. 

8.  Comments:  This  column  is  provided  for  a short 

term  description  of  the  module. 
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B.2  SOFTWARE  MODULES  LISTED  BY 
REFERENCE  NUMBER 

This  table  includes  all  of  the  software  modules  listed 
in  numerical  order  by  module  reference  number. 

100  SERIES:  Basic  transfer  and  input/output 

operations 

200  SERIES:  Basic  mathematical  operations 

300  SERIES:  Initialization  and  output  configuration 

routines 

400  SERIES:  Radar  block  simulation  modules  (black 

box  models) 

500  SERIES:  Radar  system  modules  (radar  system 

models) 


Table  B.2-la  100  SERIES  MODULES 


Parameter  Initialization 
Parameter  Initialization 
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Table  B.2-lb  200  SERIES  MODULES 


201 

DFT 

231 

RSHFTS 

202 

ZFFT 

232 

RSHFTS 

1 

203 

ZIFFT 

233 

DFTRF 

204 

CONV 

234 

DFTF0 

1 

205 

CONVMP 

235 

ADDRND 

206 

DIVA 

236 

ADDRND 

207 

ADDA 

237 

ADRNDC 

\ 

208 

CUMDIS 

239 

ADDA 

209  CUMDIS 

210  OUTCUM 

211  OUTCUM 

212  PDF 

213  PDF 

2 14  RNDARY 

215  RNDARY 

216  ATOD 

217  ATOD 

218  DTOA 

2 19  DTOA 

220  DTOA 

221  WEITRE 

222  WEITCP 

223  WEITMP 

224  SHIFT 

225  SHIFT 

226  SHIFTS 

227  SHIFTS 

228  DTOA 

229  RSHIFT 

230  RSHIFT 
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Table  B. 2-lc  300  SERIES  MODULES 

301  (SYSTEM) 

302  CLINT 

303  PTLIST 

304  PTLIST 

305  PTLIST 

306  PTLIST 

307  PLOTTR 

308  PLOTTR 

309  PLOTTR 

310  PLOTTR 

311  (BISTATIC  ANTENNA  INITIALIZATION) 

312  SPCAVG 

313  SCANNR 

I I 
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Table  B.2-ld  400  SERIES  MODULES 


401 

NONLIN 

438 

HLIM 

402 

NONLIN 

439 

HLIM 

403 

CDIGFL 

440 

DCFAR 

404 

CDFNCL 

441 

DCFAR 

405 

RDIGFL 

442 

IHLIM 

406 

RDFNCL 

443 

IHLIM 

00 

o 

<r 

FILT 

445 

IHWDET 

409 

INGTOR 

446 

IHWDET 

410 

INGTOR 

447 

IFWDET 

413 

ANTARY 

448 

IFWDET 

414 

HWDET 

449 

I SQDET 

415 

HWDET 

450 

I SQDET 

416 

FWDET 

451 

RECF 

417 

FWDET 

452 

RECFTF 

418 

SQDET 

453 

CGENSF 

419 

SQDET 

454 

PXFRM 

420 

FGENXY 

455 

CGENCW 

421 

FGENMP 

456 

LAMP RE 

422 

DIGTFL 

457 

LAMPRE 

423 

DIGTFL 

458 

LAMPCP 

425 

CGEN 

459 

CFAR 

426 

TSARY 

460 

CFAR 

427 

TSARY1 

461 

DIGFIL 

430 

MTIFLT 

462 

DIGFNC 

431 

MTIFLT 

463 

DIGFSF 

432 

MTINCL 

464 

DIPOLE 

433 

MTINCL 

465 

DIPOLE 

434 

SWPINT 

435 

SWPINT 

436 

NCSWPI 

437 

NCSWPI 

] 
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Table  B.2-le  500  SERIES  MODULES 


501  TARGET 

502  TCTNCL 

503  CLUTTR 

504  ANT PAT 

505  TSRPAT 

506  PHENC 

507  PHENC 

508  PHDEC 

509  PHDEC 

510  WVGUID 

511  IONOS 

512  ECM 

513  SCYL 

514  CHAFF 

515  CHAFF 
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APPENDIX  C 

VARIABLES  USED  IN  THE 
RADAR  SIMULATION  MODEL 


This  list  includes  all  the  system  and  module  parameters 
used  in  the  modules  in  Appendix  B.  Column  headings  in  the 
list  are  defined  as  follows: 

1.  NAME:  Name  of  variable  used  in  the  system 

model  program. 

2.  SYMBOL:  Symbol  of  variable  used  in  text. 

3.  UNITS:  Dimensions  of  variable  used  in  module 

program . 

4.  LOC:  Location  of  variable  in  module  program. 

5.  MODULE  NUMBER  and  NAME:  Modules  that  use 

variable . 

6.  NAME  IN  MODULE:  Name  of  variable  used  in  the 

module  program. 

7.  COMMENTS:  Short  description,  restrictions,  and 

default  value. 

Note  that  the  variable  names  in  column  1 are  the  sys- 
tem variable  names,  which  are  referred  to  in  Volume  1.  The 
variable  names  in  column  6 are  the  module  or  subroutine 
variable  names,  which  are  referred  to  in  Volume  2.  Table 
C-2  lists  in  alphabetical  order  those  module  names  which 
are  different  from  their  corresponding  system  names. 
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Nrt*E  MODULE  MODULE  NAME  COMMENTS 
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\£)  i— I i—1 

sf  <j-  in  m 


w 

W 

w w 

fcd 

hJ  hJ 

0 

O 

0 c 
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-SYS) 


Table  C-2  MODULE  VARIABLE  NAME/SYSTEM  VARIABLE 
NAME  CROSS  REFERENCE  TABLE 


MODULE  VARIABLE  NAME 

ANTANC 

ANTP 

ARYDEL 

AZ000 

BS  IT 

CODE 

DELTHE 

DIN 

EL000 

FB 

FBX 

FBXY 

FEXT 

FF 

FFX 

FFXY 

FSTART 

F0 

F0XMT 

IC 

ICPLXF 

ICPLXI 

IDMP 

I DMY 

IDMY1 

INORM 

MAD1 

NBITS 

NIXF 

NP 

NPT 

NPTS 

NPULS 

NSUBP 

NTYPE 

N2POW 

RN000 

RTGT0 

R000 

STARPT 

STOPT 

TCELL 


SYSTEM  VARIABLE  NAME 

AZANG  or  ELANG 
ANTAZ  or  ANTEL 
TDLINE 
AZ0 

AZBST  or  ELBST 

ICODE 

ADELT 

DFTIN 

EL0 

FBCOEF 

FBR 

FBI 

SIMBW 

FFCOEF 

FFR 

FFI 

FSTRT 

RFF0  or  SIMF0 

FSTRT 

IFCODE 

ICFOR 

ICINV 

IDDUMP 

IDDUMP 

IDDUMP 

NORMFT 

IADDI 

NBITDF 

NHIST 

NSEC 

NPTAZ,  NPTEL  or  NIMP 

NRNDPT  or  NTFN 

NIMP 

NCGPLS 

NTYPER 

N2 

RN0 

RTGT 

RSIM 

ASTRT 

ASTOP 

RNGCEL 


C-22 


p 1 


Table  C-2  MODULE  VARIABLE  NAME/SYSTEM  VARIABLE 
NAME  CROSS  REFERENCE  TABLE 
(Page  2 of  2) 

MODULE  VARIABLE  NAME  SYSTEM  VARIABLE  NAME 


THT 

TNOM 

T 0 

UEXT 

UUEXT 

WT 

XIVRNG 
XI  VST 


SHPHAS 
TCGNOM 
SHT0 
XUEXT 
UEXT 
WTX 
RNG 
ST 


ADDENDA 


A 

ACYL 

AZCHFF 

CHFAZ 

DIA 

DIPDIA 

H 

HCYL 

HCHAFF 

CHFHGT 

RCCENT 

CHFRNG 

RHO 

RHOD 

THETA 

THETA D 

TSHFT 

TSHFTD 

XL 

DIPLNG 

1 


I 


I APPENDIX  D 

TABLE  OF 

TIME  SHARING  FILES 


The  table  contained  herein  is  a list  of  the  files  in 
the  radar  simulation  model.  In  the  left  column,  the  file 
or  /catalog/file  designation  that  includes  the  password 
is  listed.  The  files  beginning  with  the  letter  "SM  are  the 
Fortran-IV  decks,  and  those  beginning  with  the  letter  "0"  mm 
are  the  object  decks.  In  the  right  column  is  the  list  of 
modules  that  are  located  in  the  corresponding  file  in  the 
left  columns. 

The  Job  Definition  Files  are  those  used  to  run  the  sample 
problems  in  Section  8 of  Volume  I and  Section  4 of  Volume  III. 


Table  D-l  TIME  SHARING  FILES 


BECAGDCU  FILES 


1 


f 


♦ LODF'EX'fRJH 
♦OLGDRtRJH 

♦SUPRTlfRJH 
*DSUP1$R JH 

♦ SCI.  I NTfF'JH 
*OCL INTfRJH 

♦■5XMTR1JRJH 

♦ OXMTRlf  F'JH 

♦ 0 ft  R 5 1 N $ R J H 
*SIMEX»RJH 
♦OSIMEXtPJH 


♦SZFFTtRJH 
*OZFFT $RJH 


*SDFI NfRJH 
* 0 B F I N I R J H 


* S X ft  T R E F R J H 
*OXMTRE$RJH 

*SPLOTSfRJH 

* 0 PLOTS# R J H 

* S II I S C 1 f R J H 
*0MI SC1FRJH 

*STSfiRl$RJH 

* 0 T 5 ft  R 1 f R J H 

* / E:  L 0 C K 1 $ R J H / S M I S C 2 1 R J H 

* / E:  L 0 C K 1 1 R J H / 0 11 1 S C £ $ R J H 

* / 6 L 0 C 1C  1 $ R J H / 5 LI  P R T E t R J H 
♦V&LOCKlFRJH/OSUPEfRJH 


lift!  I,  i 

rift  ini 

I P ft  C:  K ; F'Rft  N D ; D E:  L K X ; 6 L 0 C K D ft  T ft 
I F'ftC:K;  R R ft N D , D B L K X , B L 0 C K DATA 

CLINT 

CLINT 


F'HENC 

PHENC 

ARSI  N 
MAINE' 

MAINE 

ZFFT 

ZFFT 

D C F A R ; OlGTFLi  O'DIGFL;  INGTOR,  ft  7 0 L1  , R D I G F L ; 
DIGFIL;  ECM,  DIGFSF 

C'CFRRi  DIGTFL;  CD  I GFL;  INGTOR,  ftTGD ; RDIGFL, 
DIGFIL, ECM; DIGFSF 

FGENXV;  F'NDEC;  ABORT 
FGENXV; PHD EC, ABORT 

PTLISTi  F'LOTTR 
PTL.1ST,  F'LOTTR 

FILT, WEITRE;  CONV;  SHIFT,  C F ft  R ; LftMPCP,  ERG VCR 
FILT, WEITRE, CONV, SHIFT, CFAR; LftMPCP, ERGVCP 

RECF, PXFRM, CGENSF 
RECFj PXFRM; CGENSF 

NONLIN, HWDET;  ANT  I NT; DFT 
NONL IN; HWDET; ft  N T I NT , DFT 

RNDftRV,  CUMDIS,  RTOF'DB;  REGAIN;  ELGA  IN 
R'NDftR V;  CUMDIS;  RTOF'DB;  REGAIN;  E.LGflIN 


I 


i 


D-2 


f! 


* , 'BL  OCK 1 f R JH/S AN  T 1 1R  J H 

* / B L 0 C. K IfRJH/OANTlf R J H 

* BLOC  I 1 fR JH/S ANT 2f  RJH 
P LOCI  1 fPJH/OANT  2f  RJH 

*/BLOCf  lfRJH/STGCL  fRJH 
‘ B L C<  C K 1 $ R J h , ' 0 7 G C-  L f R J H 

*SSTSAR$RJH 

*/BL  0 C i 1 IP  JH.'OT  SARf  R JH 

* . 1 E-  L 0 CP  1 f R J H , ' 5 D I G f R J H 

* BL  OCt  if RJH  ODIGfRJH 


ANT  ARY;  F'LTFMT 
ANTARY;  F'LTFMT 

SCANNR;  SF'C  ft VG 
SCANNR,  SPC  FlVG 

TARuET  . CL  LIT  7 R.  WVGUID,  I 0 N 0 S 
TftRGETi  CLUTTR,  WVGUID;  I OLIOS 

TSRF'AT,  C GEN)  TSARV 
TSRF'AT,  COEN , TSARV 

MTIFL  1,  SWF'INT 
MTIFLT;  SWF  I NT 


* T ft  F E I N 
♦CKTAPE 
*TAFTBL 


ASSEMBLY  LONG  FROG  TO  READ  TRF'E  INTO  BUFF 
PROGRAM  TO  VERIFY  RCS  MAG  TAPES 
TAPES  SENT  TO  RADC 


JOB  DEFINITON  FILES 


RANJB 

I , <3  2 

W G D ISP 

I,  8 18 

END  1ST 

I , 6 2 

IMAGEM 

1,8  IF 

DFTPRB 

1,8  4 

CHAFFM 

I,  8.  18 

NONLNF' 

1,8  5 

BISTAT 

I,  8.  19 

BARKER 

1,8  £ 

TXCODE 

I,  8.  20 

WBTGT1/2 

1,8? 

DECODE 

I,  8.  21 

WBTGCL 

I,  8.  8 

DHF'ROB 

1 , 8.  22 

NBTGCL 

1,89 

RAIN  P B 

I,  8.  22 

WAF'C  G 

1 , 8 10 

RCSSP 

I,  8 24 

ANTNNA1 

1,8  11 

TSAR B P 

I I I,  4 1 

TSARBK 

I,  8.  12 

CGENOT 

III, 4 2 

ECMODL 

1 . 8 12 

TSARTM 

111,4.2 

DIKFIX 

1,8.  14 

TSAR JB 

III,  4 4 

C FARM  D 

I ■ 8.  15 

FSANT 

III,  4.  5 

SELECTED  BECAVUBl  FILES 


♦ _ I MAG  R 
f 0 I M A G F 


I MHGP;  ERRF'RC: 
I WAG  Pi  ERR  PRC: 


♦ •RCSMf JR/SCHAFF 

♦ , 'PC  SMI  JR/OCHAF  F 


DIPOLE;  CHAFF,  EXP  I 
DIPOLE, CHAFF, EXPI 


* ,'RC  SKf  JR/SC  VU 
*/RCSMf JR/OCYL 


SCYL, BISTGT 
SCYLi BISTGT 


♦ ' RCSMf JR/ SX SR 

♦ - PC  SMf JR/OSXSA 


GAM, BESS 
GAM; BESS 


APPENDIX  E 
TCHEBYSCHEFF  FILTER 
DESIGN  DATA 


Given  the  Tschebyschef f approximation  of  a low-pass  filter, 
i.e. 


. _ I 2 - 

|F(»|  - l + 82Tn2(w)  , 

where  F(jw)  is  the  transfer  function, 

8 is  the  ripple  factor  (i.e., 

1 

1 - 2~T  is  the  peak- to-peak 

(1  + 8 )2 

ripple) , and 

Tn(w)  is  a Tschebyschef f polynomial  of 
order  n,  i .e  . , 

Tn(«)  = cos(n  cos  w)  ^ |cu|s  1 

= cosh(n  cosh  w)}[wl>l; 
th 

then  it  can  be  shown  that  the  k pole,  S^,  of  the 
transfer  function  F(jw)  is  determined  as  follows: 


sk  ~ ^k  + >k 

= -sinh(^  sinh-1  1/6) 


sin 


(2k-l)7T 

2n 


. ,1  . . -1  . ...  (2k-l)7T 

= cosh(~  sinh  1/6)  cos  * ^ — 


where  F(jw) 


Jm 

2~ T when  o>  = 1 . 

(1+6  ) 2 


2 2 

By  setting  1 + 8 Tn  (ou)  = 2 and  solving  for  o>,  we  can 

determine  a normalization  factor  for  the  poles  such  that 
o 

F (jl)  = h,  i.e.,  the  filter  response  will  be  -3  dB  with 
respect  to  the  ripple  peak  at  <*>=1.  This  normalization 
factor  is  derived  as  follows: 


E-l 


cosh(n  cosh~^(j^)  = 1/8 
u>l  = cosh  ~ cosh”*-  1/  8 


cosh  ~ sinh-^ 


If  each  pole  is  divided  by  this  factor  the  response 
will  be  -3  db  with  respect  to  the  peak  response.  It  may  be 
noted  that, if  the  order  of  the  polynomial  is  odd,  the  peak 
response  is  identical  to  the  response  at  w = 0.  Thus,  the 
above  formula  will  provide  normalization  to  the  DC  response 
of  the  filter.  However,  if  n is  even,  the  response  at  a>  = 0 
is  identical  to  that  at  the  lower  limit  of  the  ripple;  and, 
if  it  is  desired  to  normalize  the  Tchebyscheff  response  with 
respect  to  the  DC  level  of  an  even  pole  filter,  a different 
factor  must  be  used.  This  factor  is  derived  as  follows: 


Tn“(0)  = 1 for  even  n 

l_  .....  i 2 ^ 

1 

F(j0)|  - 2 2/n\ 

1 + 8ZV(0) 

1 + 8 * 

i 2 1 

1 

|F(j®3dB^  2(1  + 8 2) 

1 + S2Tn2(jo>3dB) 

solving  for  w3dB, 

2 

0 2(1  + 8^ 
Tn  (jw3dB^  ~ s ^ 

)-l  28^  + 1 

= 82 

-1 

282  + 1 

cosh(n  cosh  w3dg)  = 

82 

1 -1 

(2 82  + 1)^ 

w3dB  = cosh  n cosh 

8 

— -1 

(1  + 8 2)^ 

= cosh  n smh 

8 

Table  B-l  contains  the  pole  locations  for  Tchebyscheff 

of  order  1 through  10  for  various  ripple  factors.  The 

2 

are  normalized  for  F (l)j  = 

p<Vak>/2- 

filters 

poles 


E-2 


03948065  .98797802  -.11457731  .89126/8/.  -.17845835  .70731417  -.22487063 


Note:  Only  the  poles  In  the  second  quadrant  arc  shown.  The  conjugate  of  all  complex  poll*.  ; s implied. 


poles  in  the  second  quadrant  are  shown.  The  conjugate  of  all  complex  poles 


Note:  Only  the  poles  in  the  second  quadrant  arc  shown.  The  conjugate  of  all  complex  pul 


01912789  .99004866  -.05551131  .89313580  -.08646089  .70879659  -.10894708  .45507542  -. 12076877  . 1 5680830 


lote:  Only  the  poles  in  the  second  quadrant  are  shown.  The  conjugate  of  all  complex  poles  Is  Implied. 


APPENDIX  F 


DIGITAL  FILTER  DESIGN 


F.l  THE  BILINEAR  Z -TRANSFORMATION 


The  bilinear  Z- transformation  is  a convenient  means  of 
converting  an  S-domain  transfer  function  into  a Z-domain 
transfer  function  which  is  compatible  with  digital  filter 
design  techniques.  Of  particular  interest  in  this  appendix 
is  the  use  of  the  Z-transform  to  obtain  high-pass  or  low-pass 
digital  filter  designs  from  normalized  (i.e. , u>  , = 1)  S- 
domain  low-pass  filter  transfer  functions.  The  Dilinear  Z- 
transformation  is  defined  by  the  following  expression: 


s - y 


2 1-Z 


-1 


1+Z' 


(F-l) 


where  T is  the  discrete  sampling  interval.  This  transforma- 
tion maps  the  entire  left-half  of  the  S-plane  into  the  area 
inside  the  Z-plane  unit  circle,  and  maps  the  right-half  S- 
plane  into  the  area  outside  of  the  unit  circle.  The  S-plane 
jcu  axis  becomes  the  unit  circle  with  ju>=0  mapped  to  Z = +1 
and  jcij=+  oo  mapped  to  Z = -1.  This  results  in  a non-linear 
warping  of  the  frequency  scale  according  to  the  relation 


(F-2) 


where  u>c  is  the  S-plane  frequency  and  cu^  is  the  correspond- 
ing Z-plane  frequency.  For  convenience,  the  parameter  ft  will 

/ wdT  \ 

be  used  to  represent  tan  ( -j-  I in  the  remainder  of  this 

appendix.  Compensation  can  be  made  for  the  warping  effect 
of  the  bilinear  Z-transform  by  prewarping  the  filter  cutoff 
frequency.  This  is  accomplished  by  solving  equation  F-2 
for  T/2  and  substituting  the  resulting  expression  into  equa- 
tion F-l. 

.-1 


s ~ n 


1-Z 


1+Z 


-l 


(F-3) 


F-l 


F.1.1  Low  Pass  to  Low  Pass  (High  Pass 
to  High  Pass)  Transformation 

In  order  to  transform  a normalized  S-plane  filter  design 
(wc=l)  to  a Z-plane  representation  with  cutoff  frequency 
cjj,  equation  F-3  is  solved  for  Z after  1.0  is  substituted 

for  u>c.  The  following  equation  gives  the  relationship  between 
the  S-plane  and  Z-plane  poles  (zeros). 


_ 1/Q  + s 

Z 1/ft  - S 


F.1.2  Low  Pass  to  High  Pass  (High  Pass 
to  Low  Pass)  Transformation 


(F-4) 


In  order  to  transform  a low  pass  S-plane  transfer  func- 
tion to  its  high  pass  complement  in  the  Z-plane,  the  follow- 
ing complementary  transformation  in  the  S-plane  is  performed 
first 


S = 1/S'  (F-5) 

where  S'  is  the  new  pole  (zero)  location  in  the  S-plane. 
When  the  substitution  specified  by  equation F -5  is  made  in 
equation  F-4  the  following  transform  is  obtained 

_ 1/fl  + 1/S' 

Z " 1/fl  - 1/S' 

S + £1 

= IT^H  * <F-6) 


F.2  DETERMINING  DIGITAL  FILTER 
COEFFICIENTS 

Once  the  transfer  function  of  a filter  has  been  deter- 
mindd  in  the  Z-plane,  it  is  a simple  procedure  to  determine 
the  feedforward  and  feedback  coefficients  of  a digital  filter. 

Assume  the  Z-plane  transfer  function  is  defined  by  the 
following  equation. 

(Z-Z01)<Z-2ot)(Z-Z02)(Z-Z02> 

H(Z)  - K,K_  * — 

1 2 <Z-V<Z-V<Z-V<Z-V 


F-2 


where  the  asterisk  indicates  the  complex  conjugate.  This 
equation  can  be  rewritten  in  the  form 


; 


i 


f 


L 


H(Z)  = Kx 


(2-20i)(z-z0I) 

<z-zpi><z-V 


K, 


(Z~ZQ2 

(Z~Zp2 


) (Z~Z02 

> (Z“Zp2 


) 

) 


(F-7) 


The  second  form  of  this  equation  indicates  that  any 
physically  realizable  transfer  function  can  be  implemented 
by  cascading  quadratic  filter  sections.  Therefore,  the 
analysis  can  be  limited  to  a single  quadratic  section  with 
no  loss  in  generality.  One  section  of  equation  F-7  is 
expanded  to  give  the  following 

9 2 

- 2ZRe{zo}  + |Z  | 

G (Z)  = K — 2 ' (F-8) 

Z2  - 2ZRe|zp}+  |Zp| 

where  Re  j } represents  the  real  part  of  { f. 

Since  this  is  a quadratic  in  Z , it  should  be  possible  to 
implement  the  function  by  a double  delay  filter  as  shown  in 
Figure  F.l-1.  The  transfer  function  for  this  filter  will  now 
be  determined  and  compared  with  equation  F-8  above. 


Eo  = 


z2  E2  + B2  \ E2  + BxE2 


E2  = Ei  + A1  \ E2  + A2  J2  E2  . 


(F-9) 

(F-10) 


By  solving  equation  F-10  for  Ez  and  substituting  into  equa- 
tion F-9  , the  following  expression  is  obtained. 


+ b2z  + 1 

A^A  - Az 


( F-ll) 


If  Bi  is  factored  from  the  numerator,  the  following  expression 
results : 

F-3 


By  comparing  F-12  with  equation  F- 8 , the  values  of  the  feed- 
forward and  feedback  coefficients  shown  in  Figure  F-l  can  be 
determined  from  the  Z-plane  pole/zero  locations: 


Bi  = 


1 Kl2 

B2  = -2B1Re(ZD) 


C = v 


A1  " 


Ao  - 


2Re(Zp) 


-|zp|2. 


(F-13) 

(F-14) 

(F-15) 

(F-16) 

(F-17) 


F.3  EXAMPLES  OF  DIGITAL  FILTER  DESIGNS 

F.3-1  Low  Pass  Filter  Design 

This  example  will  demonstrate  the  design  of  a 2-pole 
Tchebyscheff  low-pass  filter  with  a 0.5  dB  pass-band  ripple, 
a 3 dB  cutoff  frequency  of  10  Hz,  and  a pass-band  gain  of  1. 
The  sample  rate  is  100  Hz.  The  S-plane  representation  (from 
the  table  in  Appendix  E)  of  the  low-pass  filter  is 

G(S)  = (S+. 5129091- j . 7224659) (S+. 5129091+j . 7224659)  * 


Substituting  the  above  parameters  into  equations  F-2  and  F-4 
the  following  results  are  obtained. 


F-5 


— 


j 


/fd  x L80\ 

ft  = tan  ( f~  i=  tan  18  = 0.3249197,  and 


ft 

Z_  = l 7 = .6475996  + j. 3315148 

ft 

Z*  = .6475996  - j. 3315148. 

P 

The  double  zero  at  infinity  is  represented  by 


Thus,  the  transfer  function  in  the  Z-plane  is 

(Z  + 1)(Z  4-  1) 

H(Z)  = (Z-. 6475996- j. 3:  15148) (Z-. 647599 6+j . 3315148)  ’ 


From  the  roots  in  the  Z-plane,  the  filter  coefficients 
are  determined  using  equations  F-13  through  F-17: 

B1  " I?  = 1 

B2  = "2(1) (-1)  = 2 

A±  = 2(. 6475996)  = 1.295199 

A2  = -( . 64759962  + .33151482)  = -.5292872. 


The  impulse  response  of  this  filter  is  shown  in  Figure  F.3.1-1. 

F.3.2  High  Pass  Filter  Design 

This  example  problem  demonstrates  the  design  of  a 2-pole 
Tchebyscheff  high- pass  filter  with  a 100  Hz  stopband,  1-dB 


F-6 


passband  ripple  and  unity  passband  gain.  The  sampling  rate 
is  500  Hz.  Using  the  table  in  Appendix  E,  the  normalized 
low-pass  filter  poles  are  determined: 


G(S)  = 


1 

(S+. 4507 682+j .7351424) (S+.4507682-j .7351424) 


a = tan  ( 375o  X 180°\  = 0.3249197 


\ / 

Using  the  low- pass  to  high-pass  2- transform, 


S + 12 
S - 12 


the  roots  of  the  denominator  are  determined  to  be 


Zp  = .5586542  - j. 4182765 


Zp  = .5586542  + j. 4182765. 


The  two  zeros  of  the  low-pass  filter  located  at  infinity  are 
transformed  to  zeros  at  Z=1  for  the  high-pass  digital  filter. 
The  transfer  function  in  the  Z-plane  is  the  following: 


H(Z)  = 


:Z  - 1)(Z  - l) 


(Z-.5586542+j .4182765) (Z-.5586542-j .4182765) 


The  digital  filter  coefficients  are  obtained  from  equationc 
F-13  through  F-17: 


B1  = 12  - 1 


E2  = -2(1)  (1)  = -2 


C - 1 


Al  - 2(. 5586542)  = 1.117308 


A2  = -.4870497  . 


The  impulse  response  of  this  filter  is  shown  in  Figure  F.3.2-1. 

F.3.3  Integrator  Design 

This  example  will  demonstrate  the  design  of  a simple 
integrator  with  a 3 dB  cutoff  of  one  tenth  the  sampling  rate 
and  no  Z-plane  zeros.  The  S-plane  representation  of  an 
integrator  is  the  following 

G<s>  ' s-TT 

Substituting  Sp  = -1  and  12  = .3249197  into  equation  F-4,  the 
Z-plane  pole  is  determined  to  be 

Z =0.5095254. 

Since  no  Z-plane  zeros  are  to  be  simulated,  the  Z-plane 
transfer  function  is 

H(Z)  = Z-. 5095254  ' 

The  filter  coefficients  are  the  following: 

BL  = 0 


C = 1 

Al  =0.5095257 


TIME  IN 

SAMPLING 

INCREMENTS 


Figure  F. 3. 3-1  IMPULSE  RESPONSE  - DIGITAL  INTEGRATOR 


APPENDIX  G 


COMPUTER  PROGRAM 
USED  IN  READING 
RADAR  CROSS  SECTION  TAPES 


The  documentation  and  listing  of  the  computer  program 
used  in  reading  radar  cross-section  tapes  is  provided  in 
this  appendix. 

SUBROUTINE  DATAIN 

1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

DATAIN  Subroutine  Not  user  ref- 

erenced 

2.  PURPOSE: 

This  subprogram  is  used  to  read  Radar  Cross  Section 
data  types  generated  by  the  General  Dynamics/Ft.  Worth 
RCS  measurement  range. 

3.  INPUT  PARAMETERS: 

None 

4.  CALLING  SEQUENCES: 

CALL  DATAIN 

Where:  Data  is  passed  via  blank  common. 

5.  RESTRICTIONS.  REQUIREMENTS . AND  MISCELLANEOUS  DATA: 

(a)  Data  read  from  the  magnetic  tape  is  placed  in  an 
array  INDAT.  In  addition,  a status  word  IERW  is 
generated  to  indicate  errors  during  read,  if  any. 
The  specification  in  Fortran  of  the  blank  common 
area  containing  the  above  follows: 


G-l 


COMMON  INDAT(309),  IDUM(3),  IERW 


IDUM(3)  allocates  storage  for  status  words  used 
by  the  subroutine . 

(b)  The  character  sets  of  the  RCS  tapes  and  the  Honey 
well  computer  differ  in  the  representation  of 
"="  and  "+"  signs.  Therefore,  this  routine  conn- 
ects all  occurrences  of  the  above  2 characters 
into  their  appropriate  representation  in  the 
Honeywell  character  set. 

(c)  The  following  GCOS  Master  Mode  Entry  (MME)  rou- 
tines are  used  in  this  subroutine: 

GEINDS  and  GEROAD.  Reference  GCOS  manual  BR43C 
for  an  explanation  of  these  routines. 

(d)  The  file  code  allocated  for  the  magnetic  tape 
is  13.  If  a different  file  code  is  desired, 
the  variable  FA  should  be  changed  accordingly. 

(e)  This  subroutine  is  structured  to  read  records 
that  are  1850  characters  in  length. 

6.  THEORY  OF  OPERATION 


Upon  entry  to  the  subroutine,  the  necessary  control 
words  for  GEINDS  are  set  up  and  GEINDS  is  called.  GEROAD 
is  called  to  halt  execution  of  the  program  until  a record 
is  read  from  the  magnetic  tape  unit.  Upon  completion  of  the 
read,  the  status  word  returned  by  GEINDS  is  processed  to 
determine  if  an  error  occurred  and  a status  word,  IERR,  is 
generated  for  use  by  the  calling  program.  The  values  of 
IERR  are  the  following: 


IERR  = 0 
= 1 
= 2 
= 3 
= 4 
= 16 


No  error  detected 
End  of  file 

Longitudinal  parity  error  detected 
Lateral  parity  error  detected 
Error  detected  but  not  one  of  the  above 
Blank  tape  on  read 


If  the  last  word  of  the  record  is  incomplete,  the  10 
is  added  to  IERR.  If  the  record  is  short  by  more  than  1 
word,  then  20  is  added  to  IERR. 


G-2 


SVMDEF 

OftTftI  N 

20 

PLOCK 

CLfiNK  COMMON 

: 0 

INDftT 

6SS 

209 

COFFER  HPEft 

40 

I F EE 

6SS 

1 

CHAP/ WORD  XMlTTftL  STATUS 

'..0 

s r ft  r w 

CSS 

2 

STATUS  RETURN  WORD 

60 

I E R R 

CSS 

1 

ERROR  RETURN  TO  CALLING  PROGRAM 

2 0 

USE 

PREVIOUS 

8 0 

[•or  ft  i n 

5ft  VE 

ENTRV  POINT 

S»0 

LO  ft 

DCW1 

100 

STft 

DCW 

110 

L Oft 

= - 2 

120 

STft 

STftTR 

120 

M M E 

GE I NOS 

140 

PTD 

150 

ZERO 

Eft, OCM 

FILE  CODE  WORD,  Oft T ft  CONTROL  WORD 

16  0 

ZERO 

STftTR 

STATUS  RETURN  WORD 

120 

HUE 

GEROftO 

100 

L [•  ft  0 

STftTR 

190 

STAG 

STftTW 

200 

ftNft 

=0020000. OU 

210 

TZE 

ERRF 

NO  ERRORS  DETECTED  CV  TAPE  CON’P 

22  0 

CMP  ft 

=0040000, OU 

220 

TZE 

EOF 

END  OF  FILE  DETECTED 

240 

C M P ft 

=0020000, DU 

250 

TZE 

DfiLRT 

DATA  ALERT  DETECTED 

260 

LOR 

= 4,  OL 

220 

T R ft 

ERRF 

UNKNOWN  ERROR  DETECTED 

260 

EOF 

LDft 

= 1.  OL 

END  OF  FILE 

2 9 0 

EXIT 

STft 

I ERR 

3 00 

TRft 

R T URN 

: 10 

DfiLRT 

LOft 

STftTR 

DATA  ALERT 

320 

ft  tj  ft 

=0002200, OU 

330 

CMPft 

=0000200, OU 

TEST  FOR  CLANK  TAPE  ON  READ 

340 

TNZ 

PCK 

350 

LC'ft 

=16, DL 

360 

TRft 

EXIT 

3 20 

PCK 

R N ft 

=0001000, OU 

360 

TNZ 

LHTF'C 

390 

LDft 

= 2,  OL 

LONGITUDINAL  PARI  TV  CHECK 

f 00 

TRH 

ERRF 

4 10 

LftTPC 

LC’ft 

= 2,  OL 

L ft T E P ft L PARI  TV  CHECK 

4 20 

ERRF 

STft 

I ERR 

430 

ft  NO 

=0202222, OL 

TEST  RECORD  LENGTH  READ 

4 40 

STG 

I PEP 

CHAR, 'WORD  REMAINING  TO  XFER 

450 

TZE 

TEST3 

460 

LLS 

21 

ft  = N U !1  C E R OF  CHAR  IN  LAST  WORD 

420 

GRS 

21 

G = N U M B E P OF  WORDS  NOT  XM I T TED 

460 

CM  PS 

= 1,  OL 

490 

TZE 

SPEC 

I f Q = 1 ; then  ZERO  - 1 

580 

T M C 

C ft  R C K 

If  Q < 1;  then  CARRY  = 0 

516 

LOG 

=20, OL 

SHORT  PEC-  MOPE  THAN  1 WORD  DEFECTIVE 

520 

T P ft 

ENOE 

510 

C P P C K 

CMPft 

= 2,  OL 

540 

LOft 

I ERR 

550 

TNC 

SPEC 

If  A < 2;  then  CARRY  = 0 

5:6 

TESTS 

CMPft 

= 2,  OL 

5 20 

TNZ 

CHTL 

If  A ^ 3;  then  ZERO  = 0 

560 

LOft 

= 0,  OL 

590 

STft 

I ERR 

6 0 0 

T P ft 

CHTL 

C 1 0 

SPEC 

LOG 

=10, OL 

SHORT  RECORD-LAST  WORD  DEFECTIVE 

I 


I 


626 

ENDE 

ASQ 

I ERR 

636 

CHTL 

LDXB 

= 0,  DU 

640 

NWORD 

CMPX0 

= 209,  DU 

650 

T2E 

RTURN 

660 

LDA 

INDAT, 0 

620 

LDX1 

>6,  DU 

600  — 



LDQ 

=0222727272790 

690 

TESTC 

NOP 

206 

CMK 

= 012,  DL 

210 

TZE 

EQUAL 

220 

CMK 

=022. DL 

230 

TZE 

PLUS 

240 

ROTAT 

ALR 

6 

250 

SBX1 

= 1,  DU 

260 

TNZ 

TESTC 

220 

STA 

INDAT, 0 

200 

ADXB 

= 1,  DU 

290 

TRA 

NWORD 

000 

EQUAL 

ERA 

=0000066, DL 

010 

TRA 

ROTAT 

820 

PLUS 

ERA 

=0000052, DL 

036 

TRA 

ROTAT 

840 

R TURN 

RETURN 

DATAIN 

850 

FA 

BCI 

1. 000013 

86  0 

STATR 

EBSS 

2 

8 26 

DCW1 

I OTD 

INDAT, 309 

BS0 

DC  M 

BSS 

1 

890 

END 

COMBINE  REC  LENG  TEST  WITH  lOCSTAT 
CHARACTER  TRANSLATE  ♦ AND  * SIGNS 
BEGINNING  OF  WORD  LOOP 

X0  IS  INDEX  REGISTER 

BEGINNING  OF  CHAR.  LOOP 
TEST  RIGHT  MOST  CHAR  FOR  = 

TEST  RIGHT  MOST  CHAR  FOR  ♦ 

END  OF  CHAR.  LOOP 

END  OF  WORD  LOOP 
CORRECT  = SIGN 

CORRECT  + SIGN 

RETURN  TO  CALLING  PROGRAM 

BUFFER  ARRAV,  NUMBER  OF  WORDS  XMITTED 


Hr.  M il  I'll  l lx  (-) 

MSI)  M"<lrl  Nii«ticr 

I sl>  »r  i al  l In at  | ...  Data  r..«| r 

Him  wluti  all  ' J In*  r Identifiers  are  repeated) 

I I •■*|»i.'ni  v (1-L’ltF,  7 -I -Band,  l-S-R.inc,  4-C-l)and,  5-X-Band) 
HI  an«l 

( I n I /at  I hi  ( Tr m a • HCV : I -W , l -MM . 3-VM,  4-MV) 

MSI)  I 

} Full  Scale  Ran^c  Sweep  In  Inches 
LSD  I 
HI. ink 

Msn  i 

| Bee  ml  S tmber  (Progressive  from  0000) 


LSD  ' 

Data  Prefix  (/) 

MSD  Block  Number 

LSD  (00  for  this  application) 

Blank 

Blank 

Sign  I 

MSD  , 


♦ qgq-J 


In*phase  component 
In  Ml  1 llvolt  s 


Repeated  1 to  128  time* 
(90  for  this  application) 


Range  Increment  per 
sample  (In  Inches)  to 
be  marked  on  each  reel 


Quadrature  component 
In  Millivolts 


LSD  I 
Sign 
MSD  , 


Amplitude  Component 
In  Millivolts 


1821 f T 11  1 

18271 
182)  i I 
1B7'*|  x;  x 'xx 
18251  XXX  X 
1 R 2 * • X’.  X lx  X 
182’  X X X X' 
1828' \l X X X 
182^' x x xl x 

18)0)  I I 

18)1  X1  X X1  X 
18)2,  X1  X X X 
1 8 3 3 j X ; X | X X 
18)4,  Xj  X | X X 
1 8 )Sj  X ; X I X X 
I8)8,x|x;x  X1 

i8);( 

1 8 )8  X X lx  X1 
I8)q; x x ' x, 
1 840 1 x X j X X 

1841  X X|X  X| 
18-2  X X jX  X 
184)  X X X X' 
184- 

1845  X X j X Xi 
184*  X X | X | X , 

1842  X XXX 

1848  X XXX] 

1849  X XXX, 

18>o  x xx lx! 


B I ank 
Blank 

Sign 
MsD  j 

Attenuation  level 

LSD  I 
MSD  | 

J Cal  phase  on  Cal  Runs  in  Degrees 

LSD  I 
Sign 
MSD  \ 


Pitch  Angle  In  Tenths  of  Degrees 


LSD  1 
Sign 
MSD  \ 


Roll  Angle  in  Tenths  of  Degrees 


Yaw,  Arlmijth,  ir  Aspect  In  Tenths  of  Degrees 


Notes 


Al  1 Recordings  <? 

800  BPI  In  BCD  Code 
J_  Fixed  Format  Bits 

Fixed  Format  Zero 

]_  Auto  or  Preset 

Bit  Changes 

X)  Data  Bit  Changes 


Longitudinal  Parity 


3/4  Inch  gap 
3/4  Inch  gap 


Next  Record 


Fig.  1 MAGNETIC  TAPE  FORMAT  FOR  SHORT  PULSE  DATA 


COMPUTER  PROGRAM  FOR 


IMAGING  COHERENT 
SHORT  PULSE  TARGET  DATA 


I 


A printed-out  listing  of  the  computer  program  for  im- 
aging coherent  short  pulse  target  data  is  provided  in  this 
appendix.  The  computer  program  documentation  together  with  a 
sample  program  is  located  in  Volume  I,  Part  3,  Section  8.17. 


F I L£  = I MftGP 

COMMON  I NDft  T ( 209  > . 1PEF,  S T ft  T W < 2 ) . I ERW 

CHftRftC  TER  INDftTl+460,  IMDfiT2<4&0,  11)0612*460,  JNDftT4*414 
ECU  I VALENCE  < I NO  FlT  <1 ),  INDftTl >,  < INDftT<  81  v,  I NDHT2), 

♦ < INDftT<16l), INDftT2), < INDfiT<241>, INDfiT4> 

COMMON/6COCK.1/NPTSF,  ITCi  fCft V,  DELTflV,  DftTft<540) 

COMMON  /BL0CK2/  THOLO,  DRftNG,  Ift,  I TH,  ICTRL 
COMMON/£!L  OCK  9/  ICTROL  , RP,  RN,  SP,  SN,  I SE T , NE ND,  THT 
C 0NMON/8LK1/BK1 < 500) 


COMMON /BLCK11/XFI < 40, 250) 
1,  1TI,  TOi,  DELI,  DDI , XRK126) 


DIMENSION  lDfiTft<540)  , 1 DftTft 1 < 540 ) , 1 DftTftGC  540) 

1IIR<2)  , HEftDR<2>  < TRftILR<2)  , I THETft<2> 

2THDESC40)  . 1FSKP<20>  ,SVSRTT<4)  , THMISS<20> 

2TSWNDOC10) 

DIMENSION  NC40), TGTIDC10) 

CHRRRCTER  TGTID 

EQUI VfiLENCE<M2.  0K1<  1)  ) 

DR  T ft  PI.  P I 0 2 , DTR,  RTD/2  14159265.  1.  5702962,  1 74522E-02,  57.  2957£ 
DIMENSION  CTft8L<50) 

DR T R CTftEL/020, 020, 020, 052, 071. 070, 067, 066,  065,  064.  062. 

* 062,  051,  050,  047,  04b,  045,  044,  042, 042, 041,  021,  020, 

* 027, 026, 025, 024, 022, 0 22. 021, Oil, 010, 007, 006, 005, 004. 

* 002, 002, OOl, OO0,  012,  060,  022,  022/ 

NftMEL 1ST /CftRDl/  TGTID, N8RTGT 

NRMEL  I ST/CRPD2/  k.ftV,  DEL  TftV,  NWNDOS 

NftMELI ST/CRRD2V  DRRNG, CTGT, THING,  THOLD. FREQ 

NfiMELl ST/CRRD4/  NFILES  , NRPS, IR0, LRST. INOISE, NTH, M2. 1 NCX 

NftMEL 1ST/CRRD5/  M2C. NT  HO,  ROELV 
NRMELIST/CRRD6/  NFSKP. IFSKP 
N ft M E L I S T / C ft R D £■' / SVSftTQ,  TSWNDQ,  IFSEP 
REftD  <5, CftRDl) 

WRITERS, 2060 ) TGTID, NBRTGT 
) FORMftT  < ' TftRGET  *' , 10R6,  ' NBRTGT  = ' . 15) 

REftD  (5, CftRD2) 

WRITE  <6, CftRD2) 

I F < NWNDOS  GT.  0 ) GO  TO  6 
WRITE  <6, 1002) 

i F OR  MR  T < ' NWNDOS  IS  INCORRECT  ' ) 

CftLL  EXIT 
5 CONTINUE 
NWNDO  = 1 
REftD  <5,  C ft R D 2 ) 

WRITE  < 6, CfiRD2 ) 

ftLMDft  = 11.  0/FREQ 

RK  = < 2*P I )/ftLMDft 

ft<  R = ftK*DRRNG 

RCTGT  = ftOCTGT 

TH0LD1  = 20  0*ftLOG10< THOLD) 

THTEST  = 10  +THINC 
DTH  = O 5 * T H I N C 
REftD  < 5 > C ft R D 4 ) 

WRITE  <6,  C ft R D 4 > 

ZNS  = INOISE 

NPTS  = LftST  - IR0  ♦ 1 

NFTSF  = NPTS/INCX 

NTH1  = NTH  + 1 

NTHN  * 2**142 


XFQ<  40,250),  XF I I < 64,250) 
I TQ, TOO, DELQ, DDC,  XRQ<120) 

, IDftTftl <540)  , 1 D ft  T ft  Q < 

, TRftILR<2)  , I THETftC 


SVSRTT<4) 


!DRTftQ<540) 
I THETft<2) 
THMISS<20) 


TGTIDC10) 


14159265, 1.  570796 j 


74522E-02, 57.  29578/ 


NBRTGT  = ' , 15) 


H-  2 


1 r I = NTHN 
i ro-uTHu 
>:  *4  T H = NTH 
NTH02»  MI  HU, '2 
>1  r 0 2 = NTH/2 


CUL< 
610 
C ?.  0 


x ur  2 


WEIGHING  FUNCTION  D £ F I N I T ION 


M T 0 2 V = N T 0 2 - 1 
DO  1?  1=1,  NTH 

:TH  = AKM*  (NT  G2X-  I ) 
c.:th  = COS < XTH ) 

W<1)=0  0 6 + 0 52*CXTH+CXTH 
1?  CONTINUE 

MR:  I TE<£.,  7676)  « 

7676  FORMATC  WEIGHTING  FUNCT 


6 FIG 


W S’  I T E ( f,  75GCO  NRTS 
00  F ORMAT  ( ' NETS  ='  , I 5 
1 ' NPTSF  = ', 


NTH,  2NS,  NPTSF 
' NTH  = ' , 15, /'  2N 


IF  (NFILES  GT.  O . AND  NFILES  LT.  20)  GO  TO  8 
WRITE  (6, 1004) 

1004  FORMAT  C'  NFILES  INCORRECT  ') 

CALL  EXIT 
8 CONTINUE 

NWDS  = S’  0 * H R P S 
NBA  = 1052 

I F ( NRPS  . GT.  O FIND. 

1 NRPS  . LE.  6 ) GO  TO  10 

WRITE  (6,2020) 

2020  FORMFlTC'  NRPS  INCORRECT  ' ) 

CALL  EXIT 
10  CONTINUE 


READ  < 5 , C ft R D 5 > 

NPTSS  = 2*  *M2C 

WRITEC6. 4001)  NPTSS,  NTHN,  NTHO,  ADELS’ 

4O01  FORMAT  •.  ' 0 NUMBER  OF  POINTS  TRANSFORMED  IN  TIME  DOMAIN 

1 NUMBER  OF  POINTS  IN  ANGLE  TRANSFORM  = 15,/  NUH 

2H0LD5  IN  RTI  = 15,  / ’ EXPECTED  MAXIMUM  AMPLITUDE  0 


0(30 
A 00 
1000 
1010 
1 ft  2 ft 
i o : 0 

1114  0 

1 0 5 0 
1 0 6 0 
10  7 0 

1 0 e 0 

1 0 A 0 
1100 
1110 


XU  - NPTSF  + N THN*KAV 
VL  = 00 

VU  = ADELS’  + NTHN*  DELTAS’ 

READ  <5, CARDS) 

WRITE  6 , 2 0 2 0 ) < IFSKP-:  J),  J = l,  NFSKP  ) 
U 20  FORMAT-:'  FILES  TO  BE  SKIPPED  ' , ,"  ', 


JMISS  = 0 
IFILE  = 1 
ISKP  = 1 
*■  *+•**- 


READ  <5, CARDS) 

SVSATTt  1 > = S V S A T Q 
TSHNI 0 1 1 ) = T SWNDQ 
II  CONTINUE 
ITHFG  = 0 

THDES(l)  = TSWNDCK  NWNDO) 

DO  40  J=2, NTH 

THDES(J)  = THDES(J-l)  + THlNC 
40  CONTINUE 


245  THF1N  = T HOES ( N TH ) 

IF  t T HF 1 H GE  260  Ci  ) 1THFG  = 1 

WRITE  <6,  7610>  NWNDO.  SVSATT(l),  THDES<1),  THDES(NTH) 
’-'10  FORMA T < 7'  ASPECT  WINDOW  NUMBERS  12.  7 ' SVSATT  = ' , FB.  2 
1 ' WINDOW  FROM  ',F6  1,  ' TO  '.F6  2.'  DEGREES  ' 

DO  45  J = 1 , NPTSF 
DO  45  JJ  = 1. NTH 
XFO  ( JJ,  J)  = 00 
XFI  < JJ,  J)  = 00 
45  CONTINUE 

ICTRL  = 1 
I OFT  = 1 
IOFASM  = 1 
1(11  SS  = e 
NN  =0 
I TH  = 0 
ISFTST  = 0 

250  IF'  IFILE  HE.  IFSKPCISKP)  ) GO  TO  270 
WRITE  <6,6000 

6006  FORMAT  <'  FILE  TEST  START  ',) 

251  CONTINUE 

CALL  OfiTAIN 
I F < I £ R W EO  0;  GOTO  251 
CALL  ERRRRC 

I F < IERI4.  £0  1)  GOTO  255 
GOTO  251 

255  WRITE  <6.  7000)  IFILE 
7000  FORMATS  BVPASSED  FILE  ',15  ) 

265  IFILE  = IFILE  + 1 

IF*.  I S K F-  GE  NFSKP  ) GO  TO  270 
ISKP  = I SKP  + 1 
GO  TO  250 

270  WRITE  (6,6007) 

6007  FORMAT  ('  FILE  TEST  END  - FILE  PROCESS  BEGIN  SUB) 

I PASS  = 1 
IOER  = 0 
IGN  = 0 
I9TST  = 0 
IK  = 0 

275  ITH  = ITH  + 1 

278  ISRT  = 1 
I HD  =2 

260  CONTINUE 

CALL  OBTAIN 
2 9 2 I S TP  = ISRT  4 89 

IF'.  ISRT  EO.  91  ) I HD  = 1 
IF'  ISRT  EO  181)  I HD  = 2 
IF'.IERW  EO  O OR  IERW  EG. 10  OR.  1 ERW.  EO.  12)  GOTO  52 
CALL  ERRPRC 

IF  ■.  IERW,  EO  1 ) GOTO  210 
IFdERIJ  GE  20)  GOTO  70 
GOTO  215 


52  IF  < I HD  EO  2 ) GO  TO  290 
I SRTE=I SRT+22 

DEC 0D£< I NDAT1, 6710)  < HE ADR < J ) , J =1 , 2 ) , 
• < I 0 A T A I ( J ) , IDATAG(J), J = 1 SRT,  ISRTE) 

6 710  FORMAT < 2A6, 2X,  22<5X, 215,  5X>  > 


x?'-rr  = i£RTE  + i 
! ; P I £ ■=  1 SHE  *24 

DEC  ODE' INDRT2, 6 720)  < 1 Dft  Tft  I ( J > . 1 Oft  TftCi  < J ' , J = ' SR  T T - I S fir  T E ) 

6 7 1 0 F 0 ’ H h T • ,2  4 ■ 5 X , 2i  5,  5 X J ) 

r s f r r = i s f r e ♦ 1 

1 E ft  T E = I S ft  T E ♦ 2 4 

DECODE'!  I MDR  T 2 . 672CO  ( I 0 ft  T ft  I < J ) , I Dft  TftCK  J ) . J=ISftTT,  ISRTE) 

1 5 ft  f T = I S ft  T E + 1 
I SRTE  =1 SRTE + 19 

DECODE:  I NOR  T4  . t7Hi)  < I DftTftl  ( J > . I Dft  Tft  DC  J ) , J = I S ft  T T , I SRTE  ) . 
* I I H • l HD).  ( T ft ' ft  1 1 ft  ( J > , J = 1 , 2 > < TftftlLR(2),  ITHETftC  I HD) 

6710  FC'ftHftT'.  19(SX,  215.  5X). 2X.  12,  2ft6,  ft5,  I?) 

CiC'  TCI  200 
2 9 £'  C Ci  N T I N Li  E 

IF  I S ft  T EC  1 ft  HD  19TST  EC!  O)  GO  TO  79 

I £ ft  T E = I S ft  T «•  2 2 

DECODE'.  IMDftTl,  6 746)  < I DftTftl  < J),  I Oft  T ftC!  C J ) , J = ] SR  T . 1 SRTE) 

S 740  F 0 ft  IT  ft  T 1 T21i  22  ( 5X.  215,  5X)) 

I S ft  T T = I S ft  T £ + 1 
I S ft  T E = I S ft  T E ♦ 2 4 

DECODEC  INDftT2.  6720)  ( I D ft  T ft  I < J ) , IDftTftQCJ),  J = 1SRTT,  ISRTE) 

I E ft T T = I SftTE  + 1 
40  I SftTE=I  SPTE+24 

50  DECODE'  I H D ft  T 2 1 6720)  < IDftTftI < J),  IDftTftOCJ).  J=JSRTT,  ISRTE) 

I S ft  T T = I S ft  T E + 1 
ISPTEslSRTE+19 

DEC  C'DE  < I NDftT4.  6 720)  CIDftTftlCJ).  IDftTftOCJ),  J«1S£T*.  ISRTE:- 

'DO  IFC  ISTF-  GE  MUDS  ) GO  TO  220 
I SET  = ISftT  ♦ 90 
I HD  = 3 
GO  TO  200 


70  CONTINUE 

N E >-  = ft  L D < 2 1 , 15,  I ft E F ) 

IFCNSft  GT.  lOCD  GOTO  72 
74  CONTINUE 

WRITE  (6,5140)  THDESCITH),  N 6 ft 

514£i  F 0 ft  M ft  T COSHOftT  RECORD  FOUND  - ftNGLE  = " , F7.  1,  5X,  ' NOR  »',I6> 
76  CftLL  DftTftl  N 

IF'IEPW  ED.  O Oft  I £ ft:  M ED  10  Oft.  I ERW.  ED  12)  GOTO  79 
CftLL  E ft  ft  ft  PC 
I F < I E ft !•!  ED  l > GO  TO  210 
I F ( I E ft  U GE  20)  GOTO  74 
GOTO  215 

72  WRITE  (6,5160) 

5 1 6 0 E 0 ft II ft T <•  SHORT  RECORD  RESUMED  TO  BE  Gftft  NOISE  ') 

I G N = I&N  + 1 
IF  IGN  GE  2)  CftLL  EXIT 
I HD  = 2 
GO  TO  260 

79  CONTINUE 

DECODE ( IMDftTl, 6012)  191,  192 
6012  FC'ftHftT-  T21 , 15,  T241,  15) 

IF  ( 191  LT  0 ft ND  192  LT  0 ) GO  TO  60 
I F ( I 9TST  ED  1 ) GO  TO  76 
WRITE  '.6,9004) 

9004  FC'RMftT  - 0XXXXXXXXXX  FIRST  RECORD  OF  FILE  E-'ftD  XXXXXXXXXX V/) 
GO  TO  76 
60  ISR'T  = 1 
ISTF'  = 90 
I HD  = 2 
I9TST  = 


1 


2580 

IGN  = 0 

,’S90 

GO  TO  52 

26  00 

C 

(.*^10 

C 

2 6 - 0 

210 

WRITE  (6,7002)  1 F I L E 

2 6 2 0 

7002 

FORMAT  ( 7 ‘ OF  ILL  ASSUMED  AT  END  OF  F I L E < ' . 12,'  ) . . 

2-242 

1APTING  WITH  RESTART  ANGLE'  ) 

26  5b 

IF  < < THFIN  - THDESCITH)  ) . LE  THINC  ) GO  TO  96 

2662 

I OF  T = 2 

2622 

IFILE  = IFILE  ♦ 1 

2690 

J9TST  = 0 

2690 

GO  TO  278 

2790 

96 

WRITE  < 6.  5240  ) IFILE 

2 710 

5 2 4 1' 

F 0 R M ft T < ' FILE',  15,  'IS  ASSUMED  COMPLETE'  ) 

2 72  0 

IOFASM  = 2 

2 7 ? 0 

GO  TO  4O0 

2740 

C 

2750 

C 

2 76P 

215 

CONTINUE 

2 7 7 0 

WF'ITE  <6,70020  IFILE,  THETA 

2760 

700  2 

FORMAT  < ' ERROR  ATTEMPTING  TO  READ  F I LE < ' - I 2- ' ) , ASI 

27?0 

IOER  = IOER  ♦ 1 

2S00 

IF  < IOER  6E  16)  CALL  EXIT 

2S10 

GO  TO  76 

2829 

C 

2820 

C 

284  0 

220 

I6N  = 0 

2850 

A6SITH  = I THETA<1)  - ITHETA <2> 

266  0 

IF  < ft  6 S < ft  B S I T H ) GT  THTEST  ) GO  TO  120 

2870 

119 

THETA  = I TH£TA( 1)710. 

2880 

GO  TO  120 

2 090 

120 

IF  ( ITH  E 0 1 ) GO  TO  119 

290O 

T 1 = ITHETA <1)710 

2910 

T 2 = IThETA<2)710. 

2920 

D 1 = ABS  < THETA  - T1  ) 

2920 

D2  = ABS < THE  I A - T2  ) 

2940 

I F < 01  LT  D2  ) GO  TO  125 

2950 

THETA  = T 2 

2960 

GO  TO  120  f 

2 97  0 

12  5 

THETA  = T 1 

2 98  0 

120 

CONTINUE 

2990 

IF  < THETA  LT  200  AND  I THFG  EO  1 ) THETA  = 

2000 

C 

2010 

ISPTST  = ISPTST  + 1 

i 020 

IF  < ISPTST  NE  2 ) GO  TO  219 

2020 

IF'  S V S ft  T T ' 1 ) EG  I 1 ft  < 1 ) 7 1 0 OR.  S VS  AT  T < 1 ) . EC!  I I ft  < 

2040 

WRITE  <6.6015)  SVSftTT(l),  IFILE 

2050 

6015 

FORMAT  < • O « ♦ * » « ATTENUATION  LEVEL  OF  TAPE  WAS  NOT 

3060 

17  ' FILE  NO  =', 15  ) 

2 07  0 

CALL  EXIT 

2 •-<  8 0 

2 19 

CONTINUE 

2 090 

217 

CONTINUE 

2100 

GO  TO  <224, 221), IOFT 

2110 

221 

I OF  T = 1 

2 1 2 0 

IF'  SVSATKl)  EC'  I I ft  < 1 ) 7 1 0 OP  S VSfi  T T < 1 ) EG!  I 1 A < 

2120 

WRITE  <6,  6015  > S V S A T T < 1 ) , IFILE 

2140 

CALL  EXIT 

1 . . 

900 

CONTI NUE 

2 It  d 

TEST  = THETA 

2170 

J = 0 

2180 

222 

J 3 J * 1 

2 190 

I F < ft  B S < T H D E S < J ) - RBS< THETA)  ) GT  DTH  ) GO  TO  : 

2 2 00 

NN  = NN  * J - ITH 

2 2 1 0 

C 

2220 

228 

CONTINUE 

.2  0 

C 

REPLACE  DATA  ST 


H-6 


AD-A031  480 


UNCLASSIFIED 


GENERAL  DYNAMICS  FORT  WORTH  TEX  CONVAIR  AEROSPACE  DIV  F/G  17/9 

ENDO  ATMOSPHERIC-EXO  ATMOSPHERIC  RADAR  MODELING.  APPENDICES*  A— ETC(U) 
JUN  76  R J HANCOCK*  F H CLEVELAND  F30602-73-C-0380 

RADC-TR-76-186-V0L-4-PT-1  NL 


...  to 

C 

50 

I T H * J 

3 2 SO 

IK  = J 

3220 

00  ro  240 

:zfi 

C 

2 2 9 0 

2 2 4 

C OMT ! HUE 

2 r o o 

2'  2 0 

!F<  AES  THOESdTH'  - A8S<  THE  TA)  ) LE  OTH  > GO  10  240 

2:10 

if  < r m e r ft  lt  thdesrith/  > 00  to  225 

2220 

I F < I ft  ft  S S EC  1 > GO  TO  250 

22  2 6 

WRITE  '{•.  2010)  T H 0 E S < I T H ) 

2 2*6 

3010 

FORMAT!  ' ASPECT  ANGLE  ’,F6  1,'  HISSING  ' ) 

2 256 

I MISS  = I MISS  * 1 

22b0 

JMISS  = 'MISS  ♦ 1 

2 320 

THMISS  JMISS)  = THC'ESC  I TH) 

3 380 

:f<  i miss  lt  5 » go  ro  331 

3 i 9 0 

WRITE  ' tf- . 2011) 

_ tuu 

9011 

FORMAT'.  TOO  M AN V ASPECT  ANGLES  MISSING  ' ) 

3 4 10 

CALL  EXIT 

34  20 

221 

CONT I NOE 

3 4 30 

I F ( THOESCITH)  GE  THFIN  > GO  TO  400 

2 4 4 0 

I TH  = I TH  4 1 

34^0 

IK.  = IK  ♦ 1 

“’460 

GO  TO  220 

, 4 ? 0 

c 

3 4 0 0 

235 

IF  < IPASS  £02)  GO  TO  228 

3 190 

GO  TO  250 

3500 

240 

IPASS  = 2 

3510 

HU  = HU  4 l 

3520 

c 

3530 

WRITE  ' b . 2O20)  NN,  I TH,  ( HE ADR < J).  2 = 1.  2) . I I A<1>.  < TRAILRC  J).  J*l.  2) 

3540 

1 TRAILRCO,  ITHETAC1) 

3 550 

3 0 2 0 

FOPMAT  ( ' OSNEEPC' . I 2.  ' ).  THDESC ' , I 2. ' > ' , 2A6.  2X.  1 2. 2AS, A5.  I 7 > 

3 5 S 0 

2 50 

CONTINUE 

3520 

IMISS  = 0 

3580 

GO  TO  <228. 255), IPASS 

3590 

c 

3 ? 0 0 

2 55 

CONTI  HUE 

3b  1 0 

IK  = IK  4 l 

3b20 

J F = 0 

3b30 

00  25b  J = 1 , N P T S < INCX 

3t'40 

JF  = JF  4 1 

3b  5 0 

JJ  = I R 0 4 J - 1 

3 b b 0 

XI  = 10  A T A I < J J ) 

3 C'  ? 0 

x o = j o a r h o ■;  J 3 ) 

3 S 6 0 

VFI ' IK  , JF ) = Xl/10000. 

3 1‘  9 0 

r.FC",  ll  , JF  > = 20/10000 

3 2 0 0 

2 5 b 

C OMT I HUE 

3 310 

IF'  THOES(ITH)  LT  THPIN  ) GO  To  225 

3 320 

c 

3 3 3 0 

c 

3 34  0 

4 00 

C ONTINUE 

.■’9  5 0 

WRITE  < d . 3 0 0 0 ) 

: r s 0 

2 o 0 0 

FORMAT  <•  FILE  PROCESS  E NO  ') 

3 3 " 0 

c 

3 2 8 0 

IF  ■ J M I £ £ £ 0 0 ) G 0 TO  4 S 4 

J 3 9 0 

Wf  I T£,.  £.,  2S0C1)  THMISS'.  J).  J»l.  JMISS) 

. 900 

- 

F £l  .5  ft  M J • / ' ' ♦♦HHHMMMfUM'  4 4’  4 4'  4=  4:  4 4:  4.  4 4 4 '/  4 4 4’  4 4*  4 4:  4 4 4 4 4 4 4 4 4 

b 1 0 

14  f 4 4 4 4 4 4 f 4 / i*  T m E F C • L 0 WING  ASPECT  ANGLES.  HEPE  M I SS  1 NG*  //*  (Fid 

3 *■ : 0 

i . / . , 4 4’  ♦ 4 4 4 4:  4 4 4 4 4 4 4 4 4 4 4 4 4 4 4.  4 4;  4'  4 4 4 4 4.  ’/  ' 4 4 4 4 4 4’  4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 

; 8 3 O 

2 4 4 4 4 '/  / > 

3 6 4 0 

4 c 4 

CONTI  NUE 

3 9 5 0 

c 

b 0 

c 

*****  APPLV  CORRECTION  FACTOR  ANO  WEIGHTING  FUNCTION 

3 3 0 

c 

5 9 S 0 

00  459  J = 1 . NPTSF 

2 99  0 

CRA  = A K R ♦ J -ACTGT 

H-7 


tttmmrnm 


[ 


5220 

1002 

FORMftT< 1H  . ' I TEST-'. 04) 

5219 

RE  TUP N 

5240 

550 

IR£ftD=IREftD-i 

5250 

HR  I T E >.  I U.  1008)  IREftD.  IERR1.  1 ERP2.  IERR2.  IERR4,  IERR5.  IERR6 

52i9 

1008 

F ORMftT  < 1H  .'END  OF  FILE  DE TE C T ED' , 1 4.  ' RECORDS  F-ROCESSED  BY 

5 279 

• 

' ERR8RC  ERR-'. 6<J3, IX)) 

5288  ~ • 

IFUERR  £0  7)  CALL  EXIT 

5290 

IREftD=0 

5100 

I E R R 1 * 0 

5119 

I ERR2-0 

5129 

1 ERR2=0 

5119 

I ERR4=0 

5248 

IERR5-0 

5250 

1 E RR6  = 0 

5280 

RETURN 

5220 

600 

1ERR2=IERR2*1 

5280 

WRITE< IW. 1004)  IREC.ISTftT 

5290 

1004 

FORHftTClH  , ' LONFCK-' . R3. 04) 

5400 

GO  TO  800 

5410 

650 

IERR2=IERR2*1 

5420 

NRI TE< IN. 1005)  IREC, ISTRT 

5420 

1005 

FORMftT ( 1H  . ' LftTPCK-' . R2, 04) 

5440 

GO  TO  800 

5450 

700 

IERR4=IERR4*1 

5460 

URITEC IN, 1806)  IREC.ISTftT 

5470 

1006 

FORMftT < 1 H , ' UNKERR-' , R3, 04) 

5480 

00  TO  800 

5490 

750 

IERR5=IERR5*1 

5000 

GOTO  800 

5510 

850 

NRI TE< IN. 1015) 

5520 

1015 

FORMftT < ' BLftNK  TftPE  ON  REftD-THI S JOB  NILE  TERMINATE' ) 

5520 

GOTO  550 

5540 

END 

APPENDIX  I 


TSAR  SIMULATION 
RUN  CATALOGUES 


Listings  of  all  the  computer  program  runs  for  the 
TSAR  Simulation  are  provided  in  this  Appendix.  A detailed 
explanation  of  these  tables  is  found  in  Volume  3 of  this 
series . 
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ANTENNA  LOBE  SCANNER 


OUTPUT  DATA 


The  output  data  of  the  Antenna  Lobe  Scanner  computer 
program  TSAR  Simulation  runs  catalogued  in  Appendix  I is 
provided  in  this  appendix.  The  relative  gain  and  direction 
of  the  antenna  pattern  side  lobes  are  given  for  various 
absolute  gains  of  the  main  lobe.  Listings  are  given  for  side 
lobes,  both  in  order  of  magnitude  and  in  order  of  position. 

Each  page  represents  one  run,  which  is  identified  in  the 
upper  right  hand  corner  of  the  page.  The  snub  number  (SNUMB) 
and  the  data  ;e . g.  1-22)  correspond  to  the  first  two  columns 
of  the  catalogues  in  Appendix  I.  The  pass  number  refers  to 
the  entry  number  for  that  snub  number  and  date.  For  example, 
"SNUMB=8902T , 1-22,  Pass  #2"  shown  on  page  J-3  refers  to  the 
second  entry  of  "8902T,  1-22"  on  page  1-2. 

Plots  illustrated  in  Volume  III  use  these  same 
identification  markings. 
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K.l  INTRODUCTION 


This  appendix  documents  the  results  of  a study  and 
investigation  directed  toward  the  design  of  an  Interactive 
Radar  Simulation  system.  The  Interactive  Radar  Simulation 
system  described  herein  will  significantly  improve  the  eff- 
iciency of  design  engineers  in  performing  radar  system  evalu- 
ations, in  testing  new  concepts,  and  in  predicting  the  effect 
of  subsystem  parameter  changes  on  total  system  performance. 
The  suggested  design  includes  features  that  will  streamline 
the  process  by  which  the  RADC  Radar  System  Simulation  Model 
(RADSIM)  computer  program  is  used. 

This  streamlining  consists  of  eliminating  the  menial, 
time-consuming  tasks  which  are  associated  with  the  present 
method  for  submission  of  batch  or  time  sharing  computer  jobs. 
An  example  of  this  is  the  typing  of  simulation  control  and 
module  parameter  cards  on  either  a keypunch  machine  or  a 
time  sharing  terminal.  At  the  output  end  of  the  simulation 
process,  the  digestion  of  output  data  by  the  user  is  improved 
by  the  use  of  a high  resolution  CRT  graphics  display.  This 
device  will  allow  the  plotting  of  output  data  directly  from 
the  host  computer.  In  addition,  all  information  required 
to  convert  a radar  system  block  diagram  into  a form  suitable 
for  simulation  model  usage  will  be  stored  locally  on  magne- 
tic tape  in  a catalog  form  complete  with  index  and  cross  ref- 
erences so  that  it  can  be  used  in  the  same  manner  as  a ref- 
erence document.  The  hardware  system  for  implementing  these 
improvements  will  be  referred  to  as  the  Dedicated  User  Inter- 
face (DUI).  The  Interactive  Radar  Simulator  system  described 
herein  consists  oi  the  Radar  System  Simulation  Model  (RADSIM) 
computer  program,  the  Honeywell  635  host  computer  (H635)  and 
the  Dedicated  User  Interface  (DUI). 


K.2  INTERACTIVE  RADAR  SIMULATOR 
BACKGROUND  AND  REQUIREMENTS 


In  this  section  the  background  and  shortcomings  of  the 
existing  RADSIM  computer  are  discussed  and  the  Dedicated 
User  Interface  (DUI)  is  described  in  functional  terms. 

K.2.1  Simulation  Model  Usage 

The  RADSIM  computer  program  in  its  present  configuration 
consists  of  two  segments  which  are  executed  as  separate  ac- 
tivities. These  segments  are  referred  to  as  the  Simulation 
Data  Initializer  and  the  Simulation  Controller/Modules. 

The  Simulation  Data  Initializer  serves  as  the  interface 
between  the  user  and  the  Simulation  Controller.  This  seg- 
mentation of  the  simulation  model  computer  program  was  nec- 
essary since  input  data  in  a format  convenient  to  the  user 
must  be  converted  to  a form  suitable  for  use  by  the  simula- 
tion modules.  The  user  defines  the  simulation  via  two  types 
of  punch  cards:  simulation  control  cards  and  module  para- 

meter data  cards.  The  term  punch  card  is  used  here  in  the 
sense  of  either  a physical  punch  card  or  of  a card  image 
if  the  job  is  submitted  through  the  CARDIN  time  sharing  sub- 
system of  the  H635. 

The  simulation  control  cards  determine  what  operations 
are  to  be  performed  in  the  simulation  activity.  This  includes 
the  following:  (1)  scheduling  of  modules  for  execution,  (2) 

data  transfer  to  and  from  temporary  disc  files,  and  (3)  modi- 
fication of  parameters  for  multiple  executions  of  a simulation 
model  configuration. 

The  module  parameter  data  cards  define  the  parameters  of 
the  simulation  modules  to  be  executed  in  the  simulation  ac- 
tivity. Each  module  requiring  input  data  has  a unique  name- 
list  which  contains  the  input  parameters.  Upon  completion 
of  the  Data  Initialization  activity,  the  Simulation  activity 
is  executed. 

At  the  present  time  the  output  data  of  a simulation  job 
•>  s written  onto  a permanent  disc  file  or  magnetic  tape,  or 
punched  on  cards.  The  user  then  picks  up  the  output  data  at 
the  central  site  of  the  host  computer  or  uses  the  time  sharing 
system  to  dump  data  stored  on  permfile  to  paper  tape  or  to 
listing  paper. 


K.2.2  Simulation  Model  Limitations 


In  subsection  K.2.1  the  use  of  the  existing  RADSIM  com- 
puter program  was  briefly  reviewed.  In  this  subsection  those 
aspects  of  the  computer  program  which  restrict  or  otherwise 
hinder  the  use  of  the  simulation  model  will  be  discussed.  In 
performing  a simulation  job  there  are  two  phases  which  consume 
the  majority  of  the  user's  time,  i.e.  the  operations  of  data 
input  and  data  retrieval.  In  other  words,  the  major  bottle- 
neck is  the  interface  between  the  user  and  rhe  host  computer. 
Most  of  the  time  spent  in  performing  these  tasks  is  not  en- 
gineering design  time,  e.g.  the  keypunching  of  control  and 
module  parameter  cards  or  alternatively,  typing  the  cards  in- 
to a time  sharing  file.  In  addition  to  the  typing  task,  the 
user  must  constantly  refer  to  tables  of  module  names,  refer- 
ence numbers,  parameters,  and  units  of  measure  during  the  pro- 
cess of  converting  a radar  system  block  diagram  into  a form 
suitable  for  input  into  the  RADSIM  computer  program. 

At  the  present  time,  the  output  of  a simulation  job  is 
written  onto  a permanent  disc  file  or  on  magnetic  tape,  or 
is  punched  on  cards.  This  aspect  of  the  process  of  accomp- 
lishing a simulation  is  a significant  factor  in  job  turn 
around  and  is  probably  the  most  expensive  portion  when  the 
manhours  expended  in  the  manual  handling  of  data  are  considered. 

In  order  to  summarize  the  limitations  described  above, 
the  following  example  of  a simulation  job  is  offered: 

1.  The  user  defines  the  block  diagram  complete  with 
the  parameters  of  a system  to  be  simulated. 

(Elapsed  time  = 2.0  Hours,  User  Time  = 2.0) 

2.  The  user  converts  the  block  diagram  into  a deck  of 
cards  suitable  for  input  to  the  data  loader  activity. 
(Elapsed  Time  = 2.0  Hours,  User  Time  = 2.0) 

3.  The  user  submits  the  job  to  the  host  computer  and 
the  job  is  executed.  (Elapsed  Time  = 1.0  Hours, 

User  Time  =0.25  Hours) 

4.  The  user  retrieves  the  output  data  in  the  form  of 
punch  cards  or  magnetic  tape,  transports  the  data 
to  an  off  line  plotter  system,  and  plots  the  data. 
(Elapsed  Time  = 2.0  Hours,  User  Time  = 2.0  Hours) 
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In  the  above  example  the  user  spent  6.25  hours  perform- 
ing a complete  simulation  job  of  which  only  2.0  hours  were 
used  in  performing  an  engineering  design  task.  Approximately 
607c  of  the  design  engineer's  time  was  consumed  by  non-engin- 
eering tasks.  In  the  following  subsection  a system  is  des- 
cribed from  a functional  viewpoint  which  will  improve  the 
productivity  of  the  radar  design  engineer  by  freeing  him  from 
the  menial  tasks  which  must  be  performed  in  order  to  use  the 
current  RADSIM  computer  program.  In  the  example  presented 
above,  if  the  data  were  dumped  to  paper  tape  via  an  ASR  33 
teletype,  the  time  required  would  be  even  greater  since  the 
ASR  33  data  transfer  rate  is  only  10  characters  per  second 
(cps).  Since  10  to  13  characters  are  required  to  represent 
one  floating  point  number  this  corresponds  to  one  point  of 
data  per  second  or  less.  A plot  of  3000  points  would  there- 
fore, require  almost  one  hour  to  punch  on  paper  tape. 

K.2.3  Dedicated  User  Interface 

A Dedicated  User  Interface  for  the  submission  of  simu- 
lation jobs  and  retrieval  of  output  data  should,  as  a minimum, 
have  the  following  capabilities: 

(1)  Storage  for  all  computer  program  data  required  by  the 
user  for  the  specifcation  of  parameters 

(2)  The  capability  of  producing  block  diagrams  of  the  simu- 
lation job 

(3)  Storage  for  the  job  definition  file  such  that  it  can 
be  used  in  the  construction  of  future  simulation  jobs 

(4)  A high  data  rate  channel  for  communication  with  the 
host  computer.  The  upper  bound  on  data  rate  will  be 
determined  by  characteristics  of  a voice  grade  direct 
dial  telephone  channel 

(5)  A high  resolution  CRT  display  capable  of  plotting  the 
output  data 

(6)  A device  capable  of  producing  a high  quality  paper  copy 
of  anything  presented  on  the  CRT  display. 

With  a system  having  these  capabilities,  the  procedure 
of  job  setup  will  simply  be  a matt-r  of  the  user  answering 
a sequence  of  questions  posed  by  the  DUI . When  the  job  def- 
inition is  complete,  the  user  will  have  a paper  copy  of  the 
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simulation  block  diagram  and  a magnetic  tape  cassette  con- 
taining the  job  definition  file  in  a form  compatible  with  the 
RADSIM  computer  program.  The  user  then  connects  the  DUI  sys- 
tem to  the  host  computer,  selects  CARDIN  TSS  subsystem,  trans- 
fers the  contents  of  the  cassette  to  the  TSS  current  file,  and 
gives  the  RUN  command.  The  DUI  is  then  disconnected  from  the 
host  computer. 

After  the  simulation  job  has  been  executed,  the  user  con- 
nects the  DUI  to  the  host  computer.  Again  r.  question/answer 
sequence  is  entered  where  the  system  asks  the  user  what  is  to 
be  plotted  and  some  parameters  of  the  plots.  The  data  gener- 
ated by  the  simulation  job  are  then  plotted  on  the  CRT  display. 
The  CRT  display  can  subsequently  be  reproduced  on  paper  if 
the  user  desires  a permanent  copy.  In  addition  to  immediate 
plotting,  the  user  has  the  option  of  transmitting  the  contents 
of  the  simulation  job  output  file  to  the  DUI  for  recording  on 
a magnetic  tape  cassette.  The  stored  data  can  subsequently 
be  plotted  off-line.  In  addition,  the  data  generated  by  the 
host  computer  can  be  used  as  input  data  for  programs  executed 
by  the  DUI  as  a stand  alone  computer  system. 


K. 3 DEDICATED  USER  INTERFACE  HARDWARE 


A functional  block  diagram  of  the  suggested  Interactive 
Radar  simulator  is  shown  in  Figure  K.3-1.  The  DUI  subsystem 
is  indicated  by  a dashed  line.  Two  hardware  designs  which 
satisfy  the  DUI  requirements  delineated  in  paragraph  K.2.3 
are  described  in  this  section.  The  two  designs  differ  pri- 
marily in  the  minicomputer  selected  for  the  communications 
buffer.  Design  #1  is  configured  around  a Hewlett-Packard 
2100A  minicomputer.  Design  #2  is  configured  around  a Digital 
Equipment  Corporation  PDP-11/05  minicomputer.  Design  #1 
offers  performance  superior  to  that  of  Design  #2  in  every 
respect,  but  does  so  at  a higher  cost.  The  details  of  each 
design  are  discussed  and  a comparison  made  between  them  in 
the  remainder  of  this  subsection. 

K.3.1  System  Design  #1 

The  block  diagram  of  the  Design  #1  system  is  presented 
in  Figure  K.3-2.  The  features  and  parameters  of  the  sub- 
systems which  make  up  this  design  are  discussed  in  the  follow- 
ing paragraphs. 

K.3.1.1  Hewlett-Packard  2100A  Minicomputer 

The  HP2100A  minicomputer  is  a user  microprogrammable 
general  purpose  machine.  A summary  of  the  HP2100A  parameters 
is  presented  in  paragraph  K.3.3  which  compares  the  HP  and  DEC 
minicomputers . 

The  memory  requirements  for  the  HP2100A  are  determined 
by  the  combined  size  (13K)  of  the  minicomputer  Basic  Control 
System  (4K) , the  DUI  executive  program  (3K),  and  the  Simu- 
lation Job  Setup  program  (6K) . Since  memory  is  supplied  in 
either  8K  or  16K  increments,  the  memory  size  chosen  for  the 
minicomputer  is  16K. 

Eleven  Input/Output  (I/O)  channel  slots  are  required  to 
connect  the  peripheral  devices  and  special  processing  modules 
to  the  HP2100A.  The  required  I/O  interface  cards  and  other 
nrocessor  modules  are  given  by  the  following  list: 

1.  12587B  EIA  RS232  Asynchronous  data  set  interface 

2.  12531C  Current  loop  interface  to  operator's 

console 
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INTERACTIVE  RADAR  SIMULATOR  SYSTEM  BLOCK  DIAGRAM 


Figure  K.3-2  DESIGN  #1  BLOCK  DIAGRAM 
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Time  base  generator  to  provide  timed 
system  interrupts 

Write  control  store  for  user  microprogram 
deve lopment 

Controller  and  interface  to  the  Versatec 
1600A  printer/plotter  (this  item  avail- 
able from  Versatec) 

General  Purpose  interface  to  hard  copy 
interface  unit 

Serial  asynchronous  interface  to 
Tektronix  4014-1  graphics  display 

General  purpose  interface  for  fabricat- 
ing interface  to  the  Remex  Cassette  sys- 
tem and  to  the  HP9865A  cassette  memory 
unit  (2  cards  required) 

General  purpose  interface  to  the  HP2748B 
paper  tape  reader 

Floating  point  hardware. 

K.3.1.2  Anderson- Jacobson  Modem 

The  A-J  Model  ADAC  1200  Data  coupler  transmits  and  re- 
ceives data  at  1200  baud  over  the  dial-up  telephone  network 
by  wires  through  a DAA  or  acoustically  via  a telephone  hand- 
set. The  ADAC  1200  is  compatible  with  the  Bell  System  202C 
Data  Set. 


'3.  12539C 

4.  12908A 

5.  C-HP21 

6.  12620A 

7.  12531C-001 

8.  12566B 

9.  12597-002 

10.  12901A 


K.3.1.3  Teletype 

Since  the  Tektronix  display  is  primarily  intended  to  be 
a simulation  data  I/O  device,  an  inexpensive  teletype  was  in- 
cluded in  the  system  to  act  as  an  operator's  console  and  also 
to  provide  a low  speed  paper  tape  punch/read  capability. 

K.3.1.4  Remex  Magnetic  Tape  Cassette 

The  Remex  RCP0B11/3  serves  as  the  bulk  storage  device 
of  the  DUI . By  choosing  a 3-deck  device,  one  deck  can  be 
permanently  allocated  for  operating  system  and  library  stor- 
age. The  remaining  two  decks  will,  in  general,  serve  the 
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same  function  as  the  paper  tape  reader  and  punch  do  in  other 
systems.  That  is,  one  will  normally  be  the  input  file  to  an 
executing  program  and  the  other  will  be  the  output  file. 

K . 3 . I . 5 H-P  Cassette  Unit 

In  order  to  allow  data  transfer  between  the  DUI  system 
and  the  RADC  HP9820A  Calculator  system,  the  following  approach 
will  be  used.  An  additional  interface  for  the  HP2100A  must 
be  fabricated  which  is  capable  of  driving  a HP9865A  cassette 
memory  unit.  In  addition,  a new  software  I/O  driver  must  be 
written  to  properly  format  the  data  written  to  the  cassette. 

Once  accomplished,  this  will  allow  the  HP9820A  Calculator  to 
read  tapes  generated  by  the  DUI  and  vice  versa. 

K.3.1.6  Tektronix  Display 

The  Tektronix  4014-1  graphics  display  is  the  primary  de- 
vice for  information  interchange  between  the  computers  and 
the  user.  This  display  is  available  with  an  extended  graphics 
mode  which  allows  4096  x 3120  points  to  be  plotted  on  the 
display  viewing  area. 

The  Tektronix  4014  graphics  display  is  based  on  a direct 
view  storage  tube  and  requires  no  local  mass  storage  for  dis- 
play refresh.  In  addition  to  the  store  mode,  the  Tektronix 
display  has  a non-store  display  capability  which  is  called 
the  write-thru  mode.  Data  presented  to  the  user  via  this 
mode  is  not  stored  on  the  tube  face  and  appears  only  as  long 
as  the  computer  commands  it  to  be  displayed.  Data  presented 
in  this  mode  must  be  refreshed  approximately  thirty  times  per 
second . 

A capability  is  provided  by  Tektronix  for  scanning  the 
storage  tube  surface  and  outputting  a video  signal  which  is 
proportional  to  intensity.  This  output  video  can  then  be 
processed  and  printed  by  the  Versatec  printer/plotter . The 
time  required  to  scan  the  display  face  is  ten  seconds. 

K.3.1.7  Hard  Copy  Device 

The  Versatec  1600A  printer/plotter  unit  will  provide  for 
the  generation  of  paper  copies  of  data  presented  on  the  Tek- 
tronix display.  The  Versatec  printer /plotter  is  capable  of 
plotting  1600  points  across  ten  inches  at  the  rate  of  120 
scans/second.  In  addition  to  its  capabilities  as  a plotter, 
the  1600A  is  useful  as  a high  speed  line  printer  capable  of  •; 
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"J 00  lines/minute.  The  type  font  is  a 16  x 16  dot  matrix. 

The  96  ASC1T  character  set  is  standard. 

The  alternative  to  specifying  the  1 600A  would  be  to  use 
the  1600  plotter  unit  and  then  perform  the  ASCII  character  to 
dot  matrix  conversion  in  the  DUI  minicomputer.  This  approach 
would  be  more  expensive  and  would  increase  the  processing  load 
on  the  DUI  minicomputer. 

The  additional  one  line  buffer  option  is  suggested  so 
that  a new  line  of  data  can  be  transmitted  to  the  printer/ 
plotter  while  the  previous  line  is  being  plotted  or  printed. 

K.3.1.8  Hard  Copy  Interface  Unit 

The  hard  copy  interface  unit  converts  the  video,  pro- 
duced by  scanning  the  graphics  display  storage  surface,  into 
a digital  form  suitable  for  input  to  the  Versatec  printer/ 
plotter.  The  operation  of  this  unit  is  described  in  the 
following  paragraphs. 

The  video  from  the  Tektronix  display  is  passed  through  a 
threshold  device  which  outputs  a binary  "1"  if  the  video  ex- 
ceeds the  threshold  and  a binary  "0"  otherwise.  The  threshold 
setting  is  adjustable  so  that  compensation  can  be  made  if 
video  noise  level  varies.  Since  the  resolution  of  the  Tek- 
tronix display  is  4096  points  and  that  of  the  Versatec  printer/ 
plotter  is  1600  points,  a reduction  in  the  resolution  of  the 
Tektronix  output  will  be  required.  Therefore,  the  output  of 
the  threshold  must  be  stretched  and  then  sampled  at  a rate 
which  corresponds  to  4096/3  = 1365  samples  per  scan  line, 
which  is  within  the  capabilities  of  the  Versatec  printer/ 
plotter.  The  sampler  output  is  accumulated  in  a 16-bit  serial 
in/parallel  out  shift  register.  When  the  shift  register  is 
filled,  the  output  word  is  passed  to  the  minicomputer  through 
a Direct  Memory  Access  (DMA)  I/O  channel. 

In  the  minicomputer,  the  incoming  data  words  from  the 
hard  copy  interface  are  stored  in  a buffer  area  which  is  1600 
bits  in  length.  After  a scan  line  is  loaded  into  the  buffer, 
the  hard  copy  interface  is  inhibited  for  the  next  two  scan 
lines.  This  must  be  done  to  maintain  the  geometric  fidelity 
of  the  Tektronix  display  in  producing  the  hard  copy,  i.e. 
both  dimensions  are  reduced  by  a factor  of  three.  During  the 
dead  time,  the  Versatec  unit  will  request  a new  scan  line 
and  the  buffer  will  be  transmitted  to  the  Versatec  at  the 
maximum  rate  of  the  DMA  channel. 
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The  following  is  a brief  analysis  of  the  data  rates  an- 
ticipated in  performing  the  data  transfer  from  the  tektronix 
graphics  display  to  the  Versatec  hard  copy  device.  The  time 
required  to  scan  the  display  face  containing  3120  lines  is 
10  seconds.  Therefore,  312  lines  per  second  are  transmitted 
to  the  hard  copy  interface.  The  number  of  samples  taken  of 
one  scan  line  is  1365,  therefore  the  16-bit  word  rate  out  of 
the  hard  copy  interface  is  given  by  the  following  calculation. 

data  rate=  (1365  samples/line) • (312  lines/sec) /(16  samples/word) 

= 26618  words/sec. 


This  data  rate  is  easily  handled  by  the  minicomputer  DMA  chan- 
nel which  is  capable  of  over  1 million  words/second.  Since 
every  third  scan  line  from  the  display  is  passed  to  the  hard 
copy  device,  the  scan  line  transmission  rate  out  of  the  hard 
copy  interface  will  be  104  scans/second.  The  Versatec  plotter 
is  capable  of  accepting  up  to  120  scans /second . 

In  producing  a hard  copy  of  the  Tektronix  graphics  dis- 
play, the  DOI  minicomputer  is  used  simply  for  buffer  stovage. 

K.3.1.9  Paper  Tape  Reader 

Since  Hewle tt-Packard  software  is  provided  on  paper  tape 
a high  speed  paper  tape  reader  is  highly  desirable.  So  an 
interface  should  be  purchased  so  that  the  HP2748B  paper  tape 
reader  currently  in  place  at  RADC  (OCSA)  can  be  used  with  the 
system. 

K.3.2  System  Design  #2 

The  block  diagram  of  the  Design  #2  system  is  shown  in 
Figure  K.3-3.  The  features  and  parameters  of  the  subsystems 
which  make  up  this  design  are  discussed  in  the  following 
paragraphs.  The  description  of  those  subsystems  which  are 
common  to  both  system  designs  will  not  be  repeated. 

Since  the  DEC  CAPS-11  operating  system  is  resident  on 
the  DEC  cassette  unit,  a high  speed  paper  tape  reader  is  not 
needed  in  this  system. 

K.3.2.1  Digital  Equipment  Corporation  PDP-11/05  Minicomputer 

The  PDP-11/05  minicomputer  is  a general  purpose  machine 
which  is  widely  used  in  scientific  applications.  A summary 
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of  the  DEC  PDP-11/05  parameters  is  presented  in  paragraph 
K.3.3  which  compares  the  DEC  and  HP  minicomputers. 


The  memory  requirements  for  the  PDP-11/05  are  determined 
by  the  combined  size  (15K)  of  the  CAPS-11  operating  system 
(4K)  the  DUI  executive  program  (4K) , and  the  Simulation  Job 
Setup  program  (7K).  The  estimated  memory  requirements  for 
the  DEC  minicomputer  are  larger  because  no  microcode  capabil- 
ity is  available.  The  memory  size  chosen  for  this  minicom- 
puter is  16K. 

Nine  Input/Output  (I/O)  channel  slots  are  required  to 
connect  the  peripheral  devices  and  special  processing  modules 
to  the  PDP-11/05.  Four  of  the  interfaces  require  a small 
peripheral  controller  (SPC)  slot.  The  remaining  five  inter- 
faces require  system  unit  (SU)  slots.  Since  the  PDP-11/95 
has  only  three  SU  slots  available,  two  I/O  extenders  are 
required  to  interface  all  the  peripherals  and  special  modules 
to  the  PDP-11/05.  The  peripheral  mounting  box  plugs  into  one 
SU  in  the  computer  mainframe  and  provides  four  SPC's.  The 
extension  mounting  box  plugs  into  the  PDP  unibus  and  provides 
an  additional  six  SU's.  The  required  I/O  interface  cards  and 
other  processor  modules  are  given  by  the  following  list: 

1.  DC11-AB  Serial  asynchronous  system  unit 

and  clock 

2.  DC11-DA  Full  duplex  serial  asynchronous  line 

unit  with  modem  control  (p^ugs  into 
DC11-AB) . 

3.  DL11-C  Current  loop  interface  to  operator's 

console . 

4.  DD11-B  Peripheral  mounting  box  (holds  4 

SPC's) . 

5.  8-bit  parallel  interface  to  TA11-AA 

DEC  cassette  (provided  with  TA11-AA 
unit) . 

6.  DR11-C  General  purpose  interface  to  HP  9865A 

cassette  unit. 

7.  BM  792-H  Bootstrap  loader. 

- cll-A  Asynchronous  interface  to  Tektronix 

4014-1  graphics  display. 


9. 

BA11-EF 

Extension  mounting  box  (holds 

6 

SU's) 

10. 

H720-E 

Power  supply  for  BA11-EF 

11. 

DR11-B 

DMA  general  purpose  interface 
hard  copy  interface  unit 

to 

the 

12. 

C-PDP-ll/DMA 

Controller  and  interface  to  the 

Versatec  1600A  printer/plotter 
(this  item  available  from  Versatec). 


K.3.2.2  TA11  DEC  Cassette  Tape  System 

The  TA11  DEC  cassette  tape  system  will  serve  as  the 
bulk  storage  device  of  the  DUI . 

K.3.3  Comparison  of  Designs 

The  important  difference  between  the  systems  are  the 
minicomputer  characteristics,  therefore  the  remainder  of  this 
paragraph  will  consider  only  the  minicomputers.  In  Table 
K.3-1,  some  parameters  and  capabilities  of  the  two  minicom- 
puters are  listed.  The  outstanding  differences  between  the 
devices  are  the  following: 

1.  The  floating  point  arithmetic  operations  of  the  DEC 
PDP  11/05  are  performed  by  subroutines  and  therefore 
are  extremely  slow  compared  to  the  HP2100A  which 
performs  these  operations  with  a combination  of 
hardware  and  microprogramming. 

2.  The  microprogram  in  the  DEC  PDP  11/05  is  not  access- 
ible to  the  user.  The  HP2100A  is  user  microprogram- 
mable . The  advantages  of  microprogramming  are: 
increased  speed,  reduction  of  memory  space  required 
for  the  program,  and  flexibility  to  adapt  the  in- 
struction set  to  a particular  application. 

3.  The  DEC  PDP  11/05,  as  configured,  can  be  programmed 
only  in  assembly  language,  whereas  the  HP2100A  can 
be  programmed  in  FORTRAN  II  and  ALGOL  as  well  as 
assembly  language.  The  ability  to  use  FORTRAN  will 
significantly  reduce  the  programming  task  and  en- 
hance stand-alone  capability. 

4.  Hewlett-Packard  maintains  a large  library  of  HP  con- 
tributed and  user  contributed  programs  and  micro- 
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Table  K.3-1  HP2100A/DEC  PDP-11/05  COMPARISON 


HP2100A 

DEC  PDP-11/05 

Maximum  Memory  (16  bits/wrd) 

32K 

28K 

Memory  Cycle  Time  (ms) 

0.980 

0.900 

Fixed  Point  Arithmetic  (ms) 

Add 

1.96 

3.7 

Multiply 

10.7 

4.3 

Divide 

16.7 

4.8 

Floating  Point  Arithmetic  (ms) 

Add 

24-60 

376 

Multiply 

33-41 

1017 

Divide 

52-56 

2169 

No  I/O  Channels 

14 

3 

User  Microprogrammable . 

Yes 

No 

Programming  Languages 

FORTRAN 

Assembly  Language 

ALGOL 

Assembly  Language 

K-16 


... 


programs  for  the  2100  series  minicomputers.  These 
programs  are  available  to  the  user  for  a small 
handling  charge. 

5.  The  HP2100A  is  more  expensive  than  the  DEC  PDP  11/05. 
The  cost  of  the  system  configured  around  the  HP 
minicomputer  is  approximately  10%  higher  than  the 
system  configured  around  the  DEC  minicomputer. 
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K.4  DEDICATED  USER  INTERFACE  SOFTWARE 


1 


The  new  software  to  be  developed  and  the  modification  to 
be  performed  on  the  existing  RADSIM  computer  program  are  des- 
cribed in  this  subsection.  Before  proceeding  with  the  soft- 
ware description,  it  is  felt  that  the  philosophy  which  guided 
the  allocation  of  tasks  between  the  DUI  minicomputer  and  the 
host  computer  should  be  discussed. 

K.4.1  Allocation  of  Processing  Tasks 

The  data  processing  task  of  the  Interactive  Radar  Simula- 
tor is  performed  by  the  DUI  minicomputer  and  the  H635  host 
computer.  The  pertinent  factors  to  consider  are: 

1.  More  storage  is  available  in  the  H635.  However,  it 
should  be  recognized  that  although  the  H635  memory 
is  quite  large,  256K,  the  percentage  of  this  avail- 
able to  a particular  user  during  normal  working 
hours  is  usually  less  than  25%. 

2.  The  Central  Processing  Unit  (CPU)  of  the  H635  is 
faster.  Although  the  H635  CPU  is  an  order  of 
magnitude  faster  than  that  of  the  DUI  minicomputer, 
the  H635  CPU  is  time  shared  by  a number  of  users. 
Therefore  the  difference  in  speed  between  the  dedi- 
cated minicomputer  and  the  time  shared  H635  will 
probably  be  less  than  a factor  of  three. 


3.  The  data  transferred  over  the  phone  line  should  be 
minimized.  Minimizing  connect  time  should  be  under- 
stood to  mean  not  only  the  problem  of  obtaining  a 
channel  into  the  time  sharing  system,  but  also  the 
latency  of  the  host  computer  in  servicing  the  TSS 
user's  requests.  Experience  has  shown  that  the  best 
approach  is  to  move  the  desired  data  out  of  the 
machine  in  the  quickest  possible  manner.  If  binary 
data  (the  most  compact  form)  is  transmitted  to  the 
DUI  this  will,  of  course,  increase  the  load  on  the 
DUI  minicomputer. 

4.  The  DUI  minicomputer  is  dedicated  to  the  performance 
of  one  function,  whereas  the  H635  is  time  shared  by 
a number  of  users. 


- 


5.  Since  the  RADC  host  computer  is  already  installed, 
its  resources  can  be  considered  to  have  a very  low 
marginal  cost  and  therefore,  in  configuring  a sys- 
tem, it  would  be  prudent  to  maximize  the  use  of 
these  resources. 

In  view  of  the  above  factors,  especially  1,  2,  and  5, 
the  portion  of  the  RADSIM  computer  program  which  performs 
the  simulation,  i.e.,  the  simulation  activity,  should  be 
executed  by  the  host  computer.  The  function  served  by  the 
Simulation  Data  Loader  will  be  moved  to  the  DUI  system.  By 
moving  the  data  initialization  function  to  the  DUI,  all  pro- 
cessing of  user  input  data  will  be  done  by  the  minicomputer 
off-line  and  the  data  in  a compact  form,  will  then  be  passed 
to  the  host  computer  over  the  phone  line.  This  will  minimize 
the  amount  of  data  passed  over  the  phone  line  during  job 
submission . 

In  the  case  of  the  plot  programs,  the  trade  offs  are  not 
as  clear-cut.  If  minimizing  connect  time  is  important,  it 
would  be  preferable  to  have  the  plotter  programs  executed  by 
the  DUI  minicomputer  with  data  from  the  host  computer  trans- 
mitted in  binary  form.  The  binary  data  would  then  be  recorded 
on  magnetic  tape  cassette  and  processed  for  plotting  by  the 
DUI  minicomputer  at  a later  time.  The  disadvantage  of  this 
approach  is  that  the  Tektronix  plotter  software  packages 
(Terminal  Control  System  (TCS),  Advanced  Graphics-II  (AG-II) 
and  Cal  Comp  Preview)  require  more  core  space  than  is  avail- 
able in  the  DUI  minicomputer.  Therefore,  a plotter  program 
executed  by  the  DUI  minicomputer  could  not  have  all  the  fea- 
tures provided  in  the  Tektronix  software  package.  An  addition- 
al advantage  of  this  approach  is  that  since  job  output  data 
are  stored  on  tape  cassette,  the  data  can  be  reused  a number 
of  times  in  any  manner  desired  by  user. 

On  the  other  extreme,  the  simulation  output  data  would 
be  processed  by  plot  programs  in  the  H635.  The  data  trans- 
mitted to  the  DUI  would  then  be  in  the  form  of  commands  to 
the  Tektronix  display.  The  disadvantage  of  this  is  that 
connect  time  to  the  host  computer  will  be  longer  and  the  data 
received  will  be  usable  by  the  Tektronix  display  only.  Rather 
than  compromise  between  these  two  approaches,  the  system 
should  be  configured  such  that  either  technique  can  be  used. 

If  the  user  wants  only  a basic  plot  capability,  the  plot  pro- 
gram in  the  DUI  minicomputer  will  suffice.  For  more  compli- 
cated work,  the  user  will  have  the  full  Tektronix  software 
package  in  the  H635  available  to  him. 
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In  addition  to  the  above,  a mode  of  operation  should  be 
included  whereby  the  user  can  have  tabular  data  transmitted 
to  the  DUI  under  format  control  and  printed  by  the  hard  copy 
device.  This  will  be  useful  for  retrieving  tabular  output 
data  or  program  listings.  If  the  user  desires  only  to  view 
the  data  and  requires  no  hard  copy,  then  the  data  can  be  dis- 
played on  the  Tektronix  display.  The  remainder  of  the  sub- 
section contains  a description  of  the  new  software  to  be 
developed  or  purchased  and  the  modification  to  be  performed 
on  the  existing  radar  system  simulation  model. 

K.4.2  Communications  Buffer  Executive 

The  communications  Buffer  Executive  program  (EXEC)  exe- 
cuted by  the  DUI  minicomputer  serves  as  a "data  traffic  cop." 
The  EXEC  supervises  all  input/output  operations,  provides 
linkage  to  job  set  up  routines  and  plot  programs,  determines 
which  tasks  have  the  highest  priority,  handles  user  requests/ 
answers  and  system  answers/requests. 

K.4.3  Radar  System  Simulation  Model 

The  modifications  to  be  performed  on  the  Radar  System 
Simulation  Model  (RADSIM)  are  the  following: 

1.  Develop  a new  Simulation  Data  Loader  computer 
program  which  will  accept  the  job  definition  file 
from  the  DUI  system  and  load  the  data  into  a 
permfile . 

2.  Modify  the  existing  Simulation  Data  Loader  computer 
program  to  write  its  output  data  into  a permfile 
instead  of  to  system  allocated  public  disc  space. 
This  will  be  done  so  that  RADSIM  can  still  be  avail- 
able to  non-DUI  users. 

3.  Develop  a new  Data  Retrieval  computer  program  which 
reads  simulation  output  data  from  permfile,  pro- 
cesses it  according  to  mode  selected  by  the  user  and 
then  transmits  the  data  to  the  DUI  system.  The 
Terminal  Control  System  and  Advanced  Graphics-II 
plot  programs  available  from  Tektronix  should  be 
included  in  this  computer  program. 

The  portion  of  RADSIM  which  performs  the  simulation  is 
unchanged.  Therefore,  it  is  not  degraded  in  any  way  by  the 
addition  of  the  DUI  system. 
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K.4.4  Plot  Programs 

The  programs  for  converting  the  simulation  output  data 
to  a form  compatible  with  the  graphics  display  are  available 
from  Tektronix.  Some  modification  of  the  Tektronix  software 
to  be  installed  on  the  H635  will  be  required  since  it  is 
written  for  computers  having  8-bit  bytes  and  the  RADC  com- 
puter utilizes  9-bit  bytes  when  operating  in  TSS  mode. 

The  plot  program  for  installation  in  tne  DUI  minicom- 
puter is  also  available  from  Tektronix.  This  program  should 
require  no  modification. 

The  program  for  generating  a hard  copy  of  the  Tektronix 
display  will  have  to  be  developed.  This  program  will  issue 
commands  to  the  graphics  display  to  start  the  read  scan  of 
the  storage  surface,  accept  output  data  from  the  hard  copy 
interface  and  store  it  in  a buffer  area,  and  output  the  data 
to  the  Versatec  printer/plotter  at  the  appropriate  time. 

Since  the  Versatec  printer/plotter  is  essentially  slaved 
to  the  Tektronix  graphics  display  for  plot  hard  copy  generation 
there  will  be  no  plotter  software  required  for  the  Versatec 
printer /plot ter . 


K.4.5  Simulation  Job  Setup 
Program 

The  computer  program  for  the  minicomputer  which  assists 
the  user  in  the  setting  up  of  a simulation  job  is,  for  the 
most  part,  composed  of  four  subprograms.  These  are  briefly 
described  in  the  following  paragraphs. 

K.4.5.1  Job  Setup  Controller 

This  subprogram  serves  as  the  supervisor  of  all  opera- 
tions performed  during  the  job  setup  procedure.  The  Job 
Setup  Controller  asks  questions  of  the  user,  accepts  and 
processes  the  user  answers,  and  passes  commands  to  the  appro- 
priate subprograms. 

K.4.5. 2 Block  Diagram  Generator 

This  subprogram  processes  the  information  provided  by 
the  user  and  outputs  the  necessary  commands  to  cause  a block 
diagram  to  be  formed  on  the  graphics  display.  The  user  pro- 
vides the  tabular  data  for  identification  of  modules  through 
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the  graphics  display  keyboard  and  the  positional  data  for  the 
location  of  blocks  through  the  graphics  display  cursors. 

K . 4 . r> . '1  .Job  Definition  File  Generator 

The  input  data  to  this  subprogram  is  basically  the  same 
as  that  provided  the  Block  Diagram  Generator,  but  the  output 
is  in  a form  that  is  compatible  with  the  simulation  program 
executed  by  the  host  computer.  The  structure  of  this  data 
will  be  the  same  as  the  output  data  from  the  current  Data 
Initialization  activity.  The  Job  Definition  File  which  is 
the  output  from  this  subprogram  is  written  onto  a magnetic 
tape  cassette.  The  tape  cassette  will  then  be  the  user’s 
record  of  the  input  data  for  the  simulation  to  be  performed. 

A pictorial  representation  of  the  recommended  data  structure 
on  the  tape  is  shown  in  Figure  K.4-1.  The  simulation  job 
represented  in  this  figure  is  composed  of  two  configurations. 
The  first  configuration  is  to  be  executed  twice  for  two  sets 
of  input  data.  The  second  configuration  is  to  be  executed 
only  once . 

K.4.5.4  RADSIM  Catalog  and  Data  Retrieval 

As  part  of  the  Job  Setup  Program  development  effort,  all 
data  associated  with  the  RADSIM  computer  program  pertinent  to 
the  process  of  job  definition  should  be  cataloged  and  stored 
on  a tape  cassette.  During  the  process  of  job  setup,  the 
user  may  at  any  time  request  that  data  from  this  file  be  pre- 
sented either  on  the  graphics  display  in  the  scratch  pad  area, 
or  printed  on  the  operator's  console. 


CFG  #1 


Figure  K.4-1  JOB  DEFINITION  FILE  STRUCTURE 


K . 5 INTERACTIVE  RADAR  SIMULATOR  OPERATION 


The  operation  of  the  Interactive  Radar  Simulator  will  be 
described  from  the  users  viewpoint  in  the  following  paragraphs. 
Four  outstanding  features  of  the  system  are  the  following: 

1.  No  punch  card  handling  required. 

2.  No  voluminous  quantity  of  reference  material 
required  to  set  up  a simulation  job. 

3.  No  trips  to  host  computer  to  retrieve  data  from 
simulation  runs. 

4.  Permanent  records  of  input  job  definition  data  and 
simulation  output  data  are  stored  in  magnetic  tape 
cassettes  which  are  compact  and  easy  to  handle. 

The  operation  of  the  system  will  be  described  in  terms 
of  the  phases  through  which  a simulation  job  passes.  The 
phases  of  operation  of  the  interactive  radar  simulation  are 
shown  in  Figure  K.5-1. 

K.5.1  Phase  I-Simulation  Job  Setup 

This  phase  is  done  off-line,  i.e.  the  DUI  system  is  not 
connected  to  the  host  computer.  The  intent  in  this  phase  is 
to  have  all  information  stored  locally  that  the  user  must 
know  in  order  to  set  up  a simulation  job.  The  storage  media 
for  the  data  will  be  the  RADSIM  catalog  cassette  tape.  The 
user  can,  at  any  time,  call  up  various  lists  of  information 
which  includes  the  following: 

1.  Modules  available. 

2.  Input  parameters  for  each  module. 

3.  Data  storage  arrays  used  by  each  module. 

4.  The  core  required. 

5.  Input  parameter  units. 

6.  The  default  values  of  parameters. 

During  the  job  setup,  the  display  area  of  the  Tektronix 
graphics  display  is  divided  into  two  areas,  a scratch  pad 
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Figure  K.5-1  INTERACTIVE  RADAR  SIMULATOR  FUNCTIONAL  BLOCK  DIAGRAM 


area  for  the  system  questions  and  user  answers,  and  an  area 
for  display  of  the  block  diagram.  The  block  diagram  is 
written  in  write  store  mode.  System  questions  and  user  an- 
swers are  written  in  wri te- through  mode.  The  proposed  display 
layout  is  shown  in  Figure  K.5-2.  In  the  figure  a partially 
completed  block  diagram  is  shown.  The  module  names  and  ref- 
erence numbers  are  placed  inside  the  blocks  and  the  input/ 
output  data  storage  arrays  are  shown  over  the  lines  drawn  to 
connect  the  modules.  On  the  right  side  of  the  figure  the 
system  has  asked  the  user  to  provide  a value  for  the  parameter, 
LSB,  which  is  associated  with  the  ATOD  module.  The  user's 
answer,  0.0020,  is  shown  below  the  system  question.  The  a- 
mount  of  data  which  can  be  simultaneously  displayed  in  the 
scratch  pad  area  is  limited  to  approximately  32  characters. 

This  restriction  is  due  to  the  need  for  refreshing  the  dis- 
played information  at  a rate  of  30  times  per  second.  The 
maximum  data  rate  into  the  graphics  display  is  960  characters/ 
second.  Therefore,  only  32  characters  can  be  refreshed  at 
the  rate  of  30  times  per  second. 

As  the  simulation  job  is  defined,  the  information  re- 
quired by  RADSIM  for  performing  the  simulation  is  written  on 
a magnetic  tape  cassette.  When  the  user  has  completed  the 
Simulation  Job  Setup  procedure,  the  complete  description  of 
the  job  will  be  on  a tape  cassette  and  can  be  reused  or  modi- 
fied at  a later  time.  The  tape  created  during  a job  setup 
phase  can  be  replayed  into  the  system  and  stopped  at  any 
point.  The  block  diagram  will  be  re-created  on  the  display 
up  to  the  point  at  which  the  procedure  is  stopped.  The  user 
may  increment  this  process  one  block  at  a time  by  granting 
permissions  of:  proceed,  proceed  with  write  through,  skip, 

substitute  or  insert  steps.  This  will  allow  an  existing  job 
tape  to  be  used  in  creating  a new  job  tape. 

K.5.2  Phase  II- Job  Submission 

In  this  phase,  the  user  connects  to  the  host  computer 
via  a telephone  line  and  has  the  Simulation  Data  Initializa- 
tion program  loaded  into  the  H635.  The  data  on  the  cassette 
will  then  be  transmitted  to  the  H635  and  then  written  on  a 
permfile  in  preparation  for  execution  of  the  simulation.  A 
simulation  job  is  run  as  a normal  batch  job.  This  is  necess- 
ary because  of  its  large  memory  requirements,  60K  to  75K. 
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Figure  K.5-2  DISPLAY  DURING  SIMULATION  JOB  SETUP 
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K.5.3  Phase  Ill-Simulation  Execution 

The  simulation  activity  is  performed  in  the  same  manner 
as  now,  except  the  output  data  will  be  written  to  permfile 
in  binary  form.  The  current  technique  is  to  output  the  data 
under  format  control  for  subsequent  card  punch  or  transmission 
to  a remote  terminal,  typically  a teletype.  This  will  improve 
the  efficiency  of  the  permfile  utilization  by  at  least  a fac- 
tor of  2,  i.e.  twice  as  much  data  can  be  stored  per  link  of 
disc  file  space. 

K.5.4  Phase  IV-Retrieval  of  Results 

In  order  to  retrieve  the  output  data  generated  by  the 
simulation  activity,  the  user  connects  to  the  H635  via  a time 
sharing  channel  and  has  the  data  retrieval  program  loaded. 

This  program  can  then  be  commanded  to  enter  one  of  four  sub- 
modes which  perform  the  following  operations: 

1.  Transmit  binary  data  to  DUI . 

2.  Transmit  output  data  to  DUI  under  format  control. 

3.  Process  output  data  through  the  Tektronix  software 
and  then  transmit  the  data  to  the  DUI. 

4.  Transmit  output  data  to  DUI  under  format  control 
for  listing  as  tabular  data. 

The  Dedicated  User  Interface,  depending  on  the  submode 
selected,  passes  the  incoming  data  to  the  appropriate  peri- 
pheral device  or  devices. 

1.  Convert  from  binary  to  formatted  data  and  then 
store  on  magnetic  tape  cassette. 

2.  Store  on  magnetic  tape  cassette. 

3.  Pass  data  to  the  Tektronix  display  for  plotting. 

4.  Pass  the  data  to  the  Tektronix  display  or  the 
Versatec  printer/plotter  for  listing. 

Submodes  one  and  two  are  provided  so  that  data  is  passed 
over  the  communication  channel  one  time  and  can  be  processed 
and  displayed  in  any  number  of  ways  by  the  DUI  system  off- 
line from  the  host  computer.  The  disadvantage  associated 
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with  this  is  that  the  DUI  plot  programs  will  not  have  as  much 
capability  as  the  H635  plot  programs. 


K.5.5  Phase  V-Autonomous  Operation 

In  addition  to  serving  as  the  user's  interface  of  the 
Interactive  Radar  Simulator,  the  DUI  system  can  be  operated 
as  a stand-alone  minicomputer  system.  RADC  personnel  who  are 
familiar  with  FORTRAN  can  easily  make  the  transition  to  the 
computer  system  configured  about  the  2100A,  i.e.  the  Design 
4 1 system.  The  transition  to  the  DEC  11/05  offered  in  Design 
#2  will  not  be  as  simple  since  the  DEC  system  will  not  support 
a FORTRAN  compiler  and  must  therefore  be  programmed  in 
assembly  language.  The  following  two  items  available  from 
Hewlett-Packard  would  significantly  enhance  the  stand-alone 
capabilities  of  the  Design  #1  system: 

HP12907A  - Fast  FORTRAN  Processor  (Scientific 

Instruction  Set).  This#  board  plugs  into 
an  I/O  slot  and  provides  13  microcoded 
subroutines  which  greatly  enhance  the 
through-put  efficiency  of  FORTRAN  and  other 
high  level  languages. 

HP12892A  - Memory  Protect  Fence.  This  board  plugs 

into  the  memory  section  and  aids  the  opera- 
ting system  in  running  multiple  programs. 

The  addition  of  a 7 or  9 track  magnetic  tape  transport 
or  a disc  subsystem  should  be  considered  if  the  number  of 
autonomous  progressing  jobs  are  sufficient  to  justify  the 
cost . 
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HEWLETT  PACKARD  9820A 
COMPUTER  PROGRAM  FOR  PLOTTING 
FILTER  TRANSFER  FUNCTIONS 

M.l  INTRODUCTION 

The  S-domain  transfer  function  for  a general  filter  can 
be  represented  by  the  ratio  of  two  polynomials  as  shown  in 
the  following  expression 

H(s)=  (S~2l)(S ~z2)(s'23)  • • • (S~ZNZ)  . SF 
(S-P1)(S-P2)(S-P3)  . . . (S-PNP) 

where : 

represents  the  ifc^  complex  zero, 

P^  represents  the  ic  complex  pole  and 
SF  is  a scale  factor. 


The  frequency  units  used  in  the  expression  above  can  be  any 
frequency;  i.e.,  KHz,  MHz,  GHz  etc.  as  long  as  all  inputs 
are  consistent  in  the  same  units. 


M. 2 PROGRAM  OPERATION 

The  filter  program  plots  the  baseband  filter  impulse 
response  based  on  normalized  pole-zero  locations.  The  pro- 
gram can  be  used  to  plot  magnitude  in  either  voltage  vs. 
frequency  or  power  in  dB  vs.  frequency.  The  theory  of  oper- 
ation is  shown  on  page  4-61  of  the  RADSIM  report  Vol.  II, 

Part  1.  The  following  are  the  user  supplied  inputs  required 
by  the  program: 

1.  BYPASS  - Flag  set  to  zero  or  one;  a value  of 

one  indicates  that  the  user  desires 
to  key  in  the  plot  parameters  such  as 
frequency  range  and  magnitude  of  im- 
pulse response  for  plotting;  a value 
of  zero  indicates  that  the  previous  set 
of  plot  parameters  are  to  be  used  for 
the  next  plot. 
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2 . FSTART 

3 . FSTOP 

4.  LOG  PLOT 

5 . TOP 

6 . BOTTOM 

7.  INPUT 

8.  NO.  POLES  - 

9.  NO.  ZEROES  - 

10.  SCALE 

11.  POLE-R 

12.  I 

13.  ZERO-R 

14.  I 


Starting  frequency  for  plot  (frequency 
Hz,  KHz,  MHz,  etc.) 

Ending  frequency  for  plot  (frequency 
Hz , KHz  , MHz  , e tc  . ) 

Flag  set  to  zero  or  one;  zero  indicates 
the  magnitude  of  the  impulse  response 
will  be  plotted  in  units  of  voltage  vs. 
frequency;  one  indicates  a plot  of  dri 
power  vs.  frequency. 

Maximum  magnitude  to  be  plotted 
(volts  or  dB) . 

Minimum  magnitude  to  be  plotted 
(volts  or  dB) . 

Flag  set  to  zero  or  one;  zero  indicates 
use  same  pole  zero  locations  as  for 
previous  plot;  one  indicates  that  new 
pole  zero  locations  are  to  be  entered. 

Number  of  poles  to  be  entered. 

Number  of  zeroes  to  be  entered. 

Scale  factor:  used  to  normalize 

magnitude  of  impulse  response  to  unity. 

Real  part  of  complex  pole. 

Imaginary  part  of  complex  pole. 

Real  part  of  complex  zero. 

Imaginary  part  of  complex  zero. 


Once  all  the  required  input  data  has  been  supplied  by 
the  user,  the  filter  impulse  response  is  plotted;  however, 
no  axes  will  be  drawn  or  labeled.  All  inputs  are  listed  on 
the  HP9820  calculator  printer  so  that  the  user  may  label  the 
plot  if  desired.  Tables  of  normalized  pole-zero  locations 
for  Chebyshev  filters  are  listed  in  Volume  IV,  Part  1, 
Appendix  E. 
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M. 3 SAMPLE  PROBLEM 


The  following  sample  problem  will  be  a plot  of  the  im- 
pulse response  of  a six  pole  Chebyshev  filter  with  0.1  dB 
ripple.  The  HP9820  calculator  printer  listing  which  shows 
user  supplied  inputs  is  shown  in  Figure  M.3-1. 

The  output  impulse  response  plotted  in  volts  vs.  fre- 
quency is  shown  in  Figure  M.3-2  and  in  dB  vs.  frequency  is 
shown  in  Figure  M.3-3.  The  axes  and  labels  have  been  done 
by  the  user  based  on  the  input  listing  shown  in  Figure  M.3-1. 


M.4  HOW  TO  USE  THE  FILTER  PROGRAM  RESULTS 
WITH  RADSIM 

The  primary  purpose  of  running  the  filter  program  is  to 
generate  accurate  input  data  for  RADSIM  and  in  particular 
for  Module  FILT  in  RADSIM.  This  subroutine  simulates  a con- 
tinuous filter  which  is  defined  by  an  S-domain  polynomial 
transfer  function. 

When  executing  FILT,  the  user  must  use  the  impulse  res- 
ponse in  terms  of  voltage  vs.  frequency.  The  amplitude  must 
be  normalized  to  unity  gain  and  the  pole-zero  locations  must 
be  modified  to  reflect  the  desired  bandwidth  of  the  filter 
to  be  used  in  the  simulation  job. 

Suppose  the  input  job  stream,  for  this  example,  has  been 
set  up  to  reflect  consistent  frequency  inputs  in  gigahertz 
(GHz)  and  time  inputs  in  nonoseconds  (ns).  The  frequency 
scale  in  Figures  M.3-2  and  M.3-3  is  now  constrained  to  units 
in  GHz.  Suppose  that  the  desired  six-pole  Chebyshev  filter 
with  0.1  dB  ripple  has  a 3 dB  bandwidth  of  50  MHz  and  unity 
gain  at  the  peak  of  the  impulse  response.  The  following 
procedure  would  be  used  to  modify  the  pole-zero  locations  to 
reflect  the  desired  50  MHz  3 dB  bandwidth  and  to  calculate  a 
scale  factor  (SF)  to  normalize  the  peak  of  the  impulse  response 
to  unity.  In  Figure  M.4-1  the  3 dB  bandwidth  extends  from 
-1  GHz  to  +1  GHz.  If  a lowpass  filter  is  to  be  simulated, 
the  3 dB  bandwidth  would  extend  from  DC  (0  frequency)  to  1 
GHz;  however,  for  a bandpass  filter,  the  bandwidth  will  be 
measured  symmetrically  about  zero.  To  modify  the  pole  lo- 
cations, the  following  procedure  is  used. 
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Figure  M.3-1  HP9820  LISTING  SHOWING  FILTER  PROGRAM  INPUTS 
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Each  normalized  pole  location  is  multiplied  by  (3  to  deter- 
mine the  new  pole  locations  which  define  the  50  MHz  filter. 

The  scale  factor,  SF,  is  defined  according  to  the 
following  expression. 


SF  = (1/A)  (/J) 


where 


A = the  amplitude  of  the  impulse  response 
P = pole  multiplication  factor 
n = number  of  poles 

The  scale  factor  for  this  example  would  be: 

SF  = (l/8.3)(.0250)6 
SF  = 2.9415  x 10'11 

If  the  filter  transfer  function  contained  zeroes,  then  the 
scale  factor  would  contain  another  factor  as  follows: 


where 


SF  = (1/A)  (>3)n//3r 


m = number  of  zeroes 


These  results  are  now  used  by  RADSIM  module  FILT. 

An  HP9820  calculator  listing  of  the  filter  program  is 
shown  in  Figure  M.4-2. 
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Figure  M.4-2  HP9820  LISTING  OF  FILTFR  PROGRAM 


MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  In  command,  control,  and  communications 
(C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  contnuni ca ti ons , electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology , 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


